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Abstract: The phase diagram of large Nf.., weakly-coupled A/" = 4 supersymmetric 
ff^ \ Yang-Mills theory on a three-sphere with non-zero chemical potentials is examined. 

In the zero coupling limit, a transition line in the /x-T plane is found, separating a 



> 

, "confined" phase in which the Polyakov loop has vanishing expectation value from a 

^ \ "deconfined" phase in which this order parameter is non-zero. For non-zero but weak 

O \ coupling, perturbative methods may be used to construct a dimensionally reduced ef- 

Q \ fective theory valid for sufficiently high temperature. If the maximal chemical potential 

exceeds a critical value, then the free energy becomes unbounded below and no genuine 
equilibrium state exists. However, the deconfined plasma phase remains metastable, 
D \ with a lifetime which grows exponentially with (not A^^)- This metastable phase 

persists with increasing chemical potential until a phase boundary, analogous to a 
spinodal decomposition line, is reached. Beyond this point, no long-lived locally stable 
^ \ quasi-equilibrium state exists. 

\ The resulting picture for the phase diagram of the weakly coupled theory is compared 

with results believed to hold in the strongly coupled limit of the theory, based on 
the AdS/CFT correspondence and the study of charged black hole thermodynamics. 
The confinement /deconfinement phase transition at weak coupling is in qualitative 
agreement with the Hawking-Page phase transition in the gravity dual of the strongly 
coupled theory. The black hole thermodynamic instability line may be the counterpart 
of the spinodal decomposition phase boundary found at weak coupling, but no black 
hole tunneling instability, analogous to the instability of the weakly coupled plasma 
phase is currently known. 
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1. Introduction 



The large N^. limit of a gauge theory is a type of classical limit 0, and greatly simpli- 
fies the structure of the theory. Generalizing QCD from three to a large number of 
colors may appear to be a drastic modification, but the Nc = oo theory exhibits many 
qualitative similarities to real QCD ^. However, there are important differences as 
well. For example, while real QCD (or any normal theory) can only develop genuine 
phase transitions in the infinite volume limit, Nc = oo theories may have phase transi- 
tions even in finite volume, because the Nc ^ oo limit acts as a thermodynamic limit. 
One well-known example of this type is the Gross- Witten transition in two-dimensional 
U{Nc) gauge theories 0]. 

More recently, large Nc gauge theories on a finite radius three-sphere at non-zero 
temperature have been studied by Sundborg and independently by Aharony et al. 

These theories were found to have phase transitions even in the limit of zero gauge 
coupling. Relevant order parameters include the Polyakov loop and the dependence of 
free energy on Nc- In the high temperature phase, the Polyakov loop expectation 
value is non-zero, the associated center symmetry is spontaneously broken, and the 
free energy scales as A*"^ as Nc ^ oo. In the low temperature phase, the Polyakov loop 
has vanishing expectation value, the center symmetry is unbroken, and the free energy 
is order one with respect to Nc {i.e., the free energy has a finite Nc ^ oo limit). In finite 
volume one cannot operationally define confinement in terms of the free energy needed 
to separate a static test quark and antiquark to infinity. Nevertheless, we will refer to 
these phases as "deconfined" and "confining," respectively, because the corresponding 
phases of SU{Nc) Yang-Mills theory in infinite volume show exactly the same features. 
The reader should bear in mind that the justification for this terminology comes only 
from the center symmetry realization and the Nc dependence of the free energy. 

In this paper we consider N' = 4: supersymmetric Yang-Mills theory (SYM) with 
gauge group SU{Nc), on a three-sphere of radius R and in the presence of non-zero 
chemical potentials associated with the global i?-symmetry. Our goal is to understand 
the phase structure of the theory as a function of both temperature and chemical 
potential. We will focus on the weakly coupled limit of the theory, but will compare 
results with those which are believed to hold in the strongly coupled limit of the theory, 
based on the study of solutions of the dual gravitational theory. 

The outline of this paper is as follows. In Section we consider the phase structure 
of A/" = 4 super- Yang- Mills theory on a sphere with i?-symmetry chemical potentials, in 
the limit of zero gauge coupling. Section ^ examines the theory with non-zero but weak 
coupling in the high temperature regime {i.e., inverse temperature small compared to 
the radius of the spatial three-sphere). In this regime, we construct a dimensionally 
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reduced effective field theory. We compute terms in the resulting effective action up to 
fourth order in the scalar fields, and evaluate the complete one-loop scalar field effective 
potential in the special case where the fields take values in the flat directions of the 
tree-level potential. Section ^ discusses the resulting high temperature thermodynamics 
and, in particular, the instability of the theory at sufficiently large chemical potential. 
We summarize the results of our weak-coupling analysis in Section ^ and compare with 
strong-coupling results obtained from the dual gravitational theory. Some possible 
extensions are discussed briefly in Section 

Appendix ^ contains details of the reduction of the zero coupling partition func- 
tion to the matrix model discussed in Section |^. Appendix |B| presents the one-loop 
diagrammatic calculations whose results are summarized in Section ^. Appendix 
contains the details of the evaluation of the one-loop scalar field effective potential, 
using a background field method, when the fields lie along flat directions. 

The remainder of this introduction summarizes previous work on the thermody- 
namics of A/" = 4 super- Yang-Mills theory on a sphere. We begin with strong-coupling 
results which emerge from the AdS/CFT correspondence (or gauge/string duality). 
This correspondence, which remains unproven but which has passed a great many 
consistency checks, originated from the congruence of two apparently different theo- 
ries: Type JIB supergravity on AdSs x and large N^. M = ^ super- Yang- Mills theory 
P, [10]. The field theory may be regarded as living on the M x S''^ boundary of AdSs 



space ||Tl|. The correspondence between bulk and boundary spaces extends to more 
general cases including ones where the boundary is 5*^ x 5*^. The field theory that 
lives on this manifold is naturally viewed as a finite temperature field theory, with the 
temperature equaling the inverse of the 5*^ circumference. 

For an 5*^ x S'^ boundary, there are two possible bulk geometries. One is (Euclidean 
signature) AdSs/Z with the discrete group Z acting freely on the AdS space, so the 
resulting manifold is smooth but topologically nontrivial. This gives rise to two spin 
structures over the manifold. The other bulk manifold with boundary x S'^ is the (Eu- 
clidean signature) Schwarzschild-AdS black hole. This manifold is topologically trivial, 
and allows a single spin structure. Thus, the sector of AdSs/Z geometry which has the 
same spin structure as the one in Schwarzschild-AdS may flop into Schwarzschild-AdS 
by a phase transition. This is the Hawking- Page phase transition and, as argued by 
Witten [0, the signature of this transition in the dual thermal field theory on x 5*^ 
should be a confinement /deconfinement phase transition. 

Various systems that generalize this idea have also been studied. These include, on 
the gravity side, the thermodynamics of rotating black holes. A rotation may be given 
either to the AdSs bulk or to the internal sphere. The former is the Kerr-AdS black 
hole and the latter, as seen from the bulk, gives rise to the Reissner-Nordstrom-AdS 
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Figure 1: Illustrations of previous results for the phase diagram of AA = 4 large N^. super- 
Yang-Mills theory as a function of temperature and chemical potential. The left hand figure 
depicts free field results from the zero coupling limit of the theory while the right hand 
figure comes from the analysis of the Einstein-Maxwell gravitational system (with three equal 
charges) which is believed to be dual to the strong coupling limit of the gauge theory. 

(RN-AdS) black hole through the usual Kaluza-Klein mechanism whereby momenta in 
compact extra dimensions appear as charges in the bulk. It is this latter system which 
is of interest here. 

The action that governs the bulk thermodynamics is the five dimensional Einstein- 
Maxwell action with negative cosmological constant (and Euclidean signature).^ This 
action has two saddle points. One is AdS space without a black hole, the other is 
the RN-AdS black hole. One may use either the Reissner-Nordstrom charge, or the 
associated chemical potential to parameterize these solutions.^ For the grand canon- 
ical ensemble, the relevant parameters are temperature and chemical potential. The 
equilibrium state of the system corresponds to whichever of the two saddle point con- 
figurations minimizes the free energy (which corresponds to the gravity action minus 

^When the rotates equally in three independent planes, the resulting thermodynamics is de- 
scribed by the Einstein- Maxwell action. For the more general case of unequal rotations, the action is 
that of five dimensional Af = 8 gauged supergravity. This is a generalization of the Einstein-Maxwell 
action and contains scalar fields which couple to the Maxwell fields. Solutions to the equations of 
motion have been obtained by Behrndt et al. ; the resulting black hole solutions are similar to the 
RN-AdS solution of the Einstein-Maxwell action. 

^Due to the analytic continuation from Minkowski to Euclidean signature, the time component 
of the Maxwell field is pure imaginary, so the parallel transporter of the gauge field around the 
boundary circle is not a pure phase, but rather a real number whose logarithm equals the chemical 
potential divided by the temperature. 
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an appropriate boundary term [|15| , p!6|). The resulting phase diagram has a transition 
hne separating these two phases, drawn as the dashed hne in the right hand diagram 
of Fig. m. This is the generahzation of the Hawking-Page phase transition to non-zero 
charge; the original Hawking-Page transition is indicated in the figure by the dot on the 
temperature axis at zero chemical potential. This transition between AdS and RN-AdS 
black hole solutions is a first order phase transition ||T^. 



Cvetic and Gubser analyzed the thermodynamic stability of the black holes |T8 
and found an instability line in the phase diagram (in addition to the Hawking-Page 
phase transition line). Beyond the instability line, the black hole that extremizes the 
Euclidean action becomes a saddle point; this corresponds to a thermodynamic insta- 
bility. The instability line is indicated by the solid line in the right hand diagram of 
Fig. |.3 

A natural question is the fate of the theory beyond this instability line. For the 
case of the RN-AdS black hole in four dimensions, Gubser and Mitra [|l^ found that the 
instability point is almost exactly where supergravity develops tachyonic modes. (There 
is 0.7% numerical discrepancy which is believed to be a numerical analysis artifact.) It 
is generally presumed that the theory enters a new phase beyond the instability line. 
But in the absence of any explicit solution corresponding to such a putative new phase, 
it is also possible that this line represents a genuine boundary to the phase diagram, 
beyond which no equilibrium state exists. This would be analogous to the situation at 
zero coupling, discussed below. In short, the fate of the theory beyond the instability 
line is not currently known. 

We now turn to previous work on the thermodynamics of weakly coupled SU (Nc) 
J\f = 4 super- Yang- Mills theory with non-zero chemical potentials for i?-charges, in the 
large Nc limit. If one considers the theory on flat space and introduces any non-zero 
chemical potential, then there is an immediate problem — no ground state (or grand 
canonical equilibrium state) exists. A chemical potential acts like a negative mass 
squared for the scalar fields. Due to conformal symmetry, this theory (on fiat space) 



•^There are three possible conserved charges corresponding to three independent rotational planes 
of S^. Fig. |l| illustrates the case with three equal charges, where the supergravity action reduces to 
Einstein-Maxwell. (Ref. uses a different parametrization in plotting the phase diagram. After 
recasting their equations for the instability line in terms of temperature and chemical potential, one 
finds a result of the form shown in Fig. |l|.) For other charge configurations, such as a single non- 
zero charge, we find a phase diagram with similar general structure, but in which the Hawking- 
Page transition line meets the instability line at a non-zero temperature (where the horizon radius 
shrinks to zero). For these more general cases, the behavior of the system at temperatures below 
that corresponding to zero horizon radius is not currently understood. We will discuss this further in 
Section |[ 
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has no mass term. Because there are flat directions in the super-potential, any non-zero 
chemical potential immediately destabilizes the theory. 

Compactifying space by replacing flat with a three-sphere of radius R avoids this 
problem. The coupling of the scalar fields to the spacetime curvature acts like a positive 
mass squared and allows one to introduce non-zero chemical potentials, provided they 
are sufficiently small. In the limit of zero gauge coupling, the maximum chemical 
potential equals the curvature-induced scalar mass, namely 1/R. For larger chemical 
potentials, the theory is again unstable and no equilibrium state exists. Hence, as 
illustrated on the left in Fig. |l], in the free field theory the line /i = is a boundary 
of the phase diagram (just like the T = line) and is not a phase transition line. 

This point of instability has been referred to as "Bose-Einstein condensation" in 
the AdS/CFT literature {e.g., Refs. |21|]) but this is highly mis-leading terminology. 
In an interacting system, Bose-Einstein condensation corresponds to spontaneous sym- 
metry breaking of a global U{1) symmetry. In a fi-T phase diagram, a Bose-Einstein 
condensation line is a phase transition line separating distinct phases — it is not a 
boundary of the phase diagram. 

In the free theory on without a chemical potential one finds 0, a "confine- 
ment /deconfinement" phase transition as the temperature is increased, as noted earlier. 
This is indicated by the dot on the temperature axis at zero chemical potential on the 
left side of Fig. |. 

Hawking and Reall considering the free theory on with a chemical potential, 
noted the phase boundary ai fi = 1 / R, but did not find any confinement /deconfinement 
transition and concluded that the free field theory in this particular system is qualita- 
tively different from the conjectured strongly coupled dual gravity system. However, 
these authors ignored the constraint of Gauss' law and the associated requirement that 
physical states be gauge invariant. It is the restriction to gauge invariant states which 
is responsible for the confinement /deconfinement transition at vanishing chemical po- 
tential [|, 0. Maint aining Gauss' law, even at zero coupling, is necessary to obtain 
valid results for the limiting weak coupling behavior of the interacting theory. In Sec- 
tion 1^ we reanalyze the free theory with non-zero chemical potentials and the Gauss 
law constraint, and show that there is a confinement /deconfinement phase transition 
line which connects the fi = 0, T ^ transition point with the zero temperature limit 
of stability point at fi = 1/R and T = 0. 

From Fig. ^ one may also note differences in the instability/phase-boundary lines. 
In the strong-coupling gravity dual, the chemical potential grows with temperature on 
the instability line and at asymptotically high temperatures, the instability line 
rises linearly. In the free field theory, the limiting chemical potential defining the phase 
diagram boundary is independent of temperature and, as we show below, is unaffected 
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by the imposition of Gauss' law. In Section |^ we consider the theory at non-zero but 
small coupling, and high temperature, and find that a thermal mass is generated in 
addition to the spatial curvature induced mass. The resulting effective scalar mass 
(squared) remains positive until one reaches a "spinodal decomposition" line on which 
the chemical potential also rises linearly at asymptotically high temperatures. This is 
discussed further in section |[ 

As this work was nearing completion, we became aware of the recent work of Basu 
and Wadia p2[ , which discusses the canonical emsenble of TV = 4 super- Yang-Mills the- 



ory with fixed i?-charges, instead of fixed chemical potentials. These authors examined 
an approximation to the zero coupling limit of the gauge invariant partition function 
with a projection onto states of a given i?-charge, and found multiple saddle points 
with differing (but always non-zero) expectation values for the Polyakov loop. They 
also constructed a simple model illustrating the effect of turning on a weak coupling, 
and discussed similarities between their model's results and properties of i?-charged 



anti-de-Sitter black holes. At first sight, the results of Ref. appear quite differ- 
ent from the picture of the A/" = 4 phase diagram which emerges from our analysis. 
However, as we discuss in Section Ol proper consideration of the relation between 



canonical and grand canonical ensembles near a first order phase transition shows that 



the qualitative results of Ref. P2] are consistent with our phase structure. 



2. Zero Coupling Limit 

Consider free A/" = 4 SU{Nc) super symmetric Yang-Mills theory on S^. This field theory 
has a global S'f/(4) /^-symmetry (denoted S'?7(4)^) and one may introduce chemical 
potentials coupled to the corresponding conserved charges. It should be emphasized 
that the "free" gauge theory considered in this section is the zero coupling limit of 
the theory while retaining Gauss' law. This is appropriate since Gauss' law holds and 
physical states are necessarily gauge invariant at all non-zero values of the coupling. In 
particular, one must have global charge neutrality of all physical states even with an 
infinitesimally small coupling constant. 

For free gauge theories on a sphere, there are two equivalent techniques for com- 
puting the resulting partition function: counting gauge invariant states directly, or 
using a suitable functional integral representation. The former method, adopted in 
Refs. 1^, |2^, 1^, takes advantage of the conformal mapping relating operators in the 



theory on to states of the theory on the sphere. One counts the number of gauge 
invariant operators of a given dimension (which maps to the energy of the state on the 
sphere), and then sums over all such operators to obtain the partition function on the 
sphere. 
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We choose to employ the functional integral approach (also discussed in Ref. 0) 
which represents the projection onto gauge invariant states via an integral over the 
time component of the gauge field. We begin, in subsection |2.1| , by reviewing the 
introduction of chemical potentials for a maximal Abelian subgroup of a non-Abelian 
global symmetry group, and then in subsection ^]2| write down the Lagrangian of the 
theory with the chemical potentials included. Using this Lagrangian, subsection |2.3| 
presents the matrix model which describes the partition function of the free theory on 
the sphere. We sketch the derivation of this matrix model, but leave details to Appendix 
0. The resulting phase diagram for the free theory is discussed in subsection and 
shown to have a phase transition line separating a "confining" low temperature or 
small chemical potential phase from a "deconfined" high temperature/large chemical 
potential phase. In the final subsection |2]^ we show that this is a first order phase 
transition in the zero coupling limit, but note that it is possible for the transition to 
become continuous at any non-zero coupling. 



2.1 Chemical Potentials for Non-Abelian Symmetries 

Consider a system with Hamiltonian H and internal symmetry group G, assumed to be 
a semi-simple compact Lie group. A Cartan subalgebra of G is also a maximal Abelian 
subalgebra, with dimension equal to the rank of the group G. Compactness implies 
that any group element may be written as the exponential of some element in the Lie 
algebra. Let U{g) be the unitary operator representing an element g of the group, and 
define 

Z{P,g) = tT[U{g)e~^^], (2.1) 

where j3 is the inverse temperature. By assumption, U{g) commutes with H for all 
group elements g. Consequently Z{(3, g) is a class function since (using trace cyclicity), 

Z{(3, g) = tr [Uir]) U{g) U{r^)-^e-^^] = tr [UimV-') e"^^] = ZiP, r^gr') , (2.2) 

for arbitrary rj in G. Any group element g is equivalent (under conjugation by group 
elements) to some element h oi a, maximal Abelian subgroup. And any element of a 
maximal Abelian subgroup may be expressed as an exponential of a sum of generators 
of a Cartan subalgebra, h = e^'^^'^T where {Qp |p=l, ■ ■ ■ ,rank(G)} are the generators 
and {7p} are real numbers. We are adopting the convention that the generators {Qp} 
are Hermitian. Therefore, Z{f5,g) may be regarded as a function of the rank(G) real 
variables {jp}, and rewritten as 

Z(/?,7p) =tr e-^V>'^'' . (2.3) 
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After an analytic continuation 7p — > —iPfip, this is the grand canonical partition func- 
tion 

/ip) = tr exp[-(3{H - /i^ Q^)] , (2.4) 

with chemical potentials {fip} associated with a maximal set of commuting conserved 
charges {Qp}- This illustrates why, given a non-Abelian symmetry group, it is natural 
to introduce chemical potentials corresponding to a Cartan subalgebra of the group. 

2.2 J\f = A Super- Yang-Mills Lagrangian with Chemical Potentials 

Four dimensional J\f = A super- Yang-Mills theory, in flat space, may be obtained from 



dimensional reduction of A/'=l super- Yang- Mills in ten dimensions [25, We con- 
sider this theory on 5*^ x S^, and include chemical potentials associated with a f/(l)^ 
maximal Abelian subgroup of the SU{4) global i?-symmetry. As mentioned in the In- 
troduction, this theory has conformal scalar-curvature coupling terms that will appear 
as mass terms for the scalar fields. 

To write the Lagrangian, we first determine how the charges Qp associated with 
the chemical potentials transform each field of the theory. The field content includes 
six scalars which we will regard as components of an antisymmetric matrix, (pij = 
1, . . . ,4) with = —(pji, four (left-handed) Weyl fermions Aj, and a vector field A^, 
all transforming in the adjoint representation of the gauge group SU{Nc). 

The scalar fields (pij are complex, and transform under the antisymmetric repre- 
sentation 6 of SU{4:)r. One may repackage these fields as a 6, defined by 

= le'^^'^ki. (2.5) 

Since the 6 representation is the complex conjugate of the 6, it is consistent to impose 
the reality condition 0*-' = More explicitly, this is 



(2.6) 



and the reality constraint leaves six (times A^^ — 1) real degrees of freedom in the scalar 
fields. To simplify later expressions, we will relabel the three independent complex 
scalar fields as 

01 = 012 , 02 = 013 , 03 = 014 ■ (2.7) 

The counting of degrees of freedom, along with transformation properties under the 
global SU{4:)ii symmetry, are summarized in Table |1|. Note that the total number 
of degrees of freedom for bosons and fermions are equal, as usual in supersymmetric 
theories. 

''See, for example, Ref. |24) for a different and perhaps more direct discussion that leads to the 
same conclusion. 
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Field D.O.F. SU{4:) 

~p 6(iV2-l) 6~ 

A, 8(iV2-l) 4 

A, 2{Nyi) 1 



Table 1: Counting of real degrees of freedom, and SU{4:)ji symmetry representations, for the fields 
oiM = A super- Yang-Mills theory. 



We choose to represent the generators of the Cartan subalgebra, in the fundamental 
representation 4, as 

= |diag(LL-L-l) , (2.8a) 

g^ = idiag(l, -1,1,-1), (2.8b) 

g^ = idiag(l, -1,-1,1). (2.8c) 

To obtain the corresponding form in the antisymmetric representation 6, one may 
consider the transformations of a rank-2 tensor under 4(8)4. Assembling the six scalars 
into a six-component vector via 

0^=(01,0t,02,0;,03,03)^, (2.9) 

the Cartan generator matrices in this representation appear as 

= diag(l, -1,0, 0,0,0) , (2.10a) 

g^ = diag(0, 0,1, -1,0,0), (2.10b) 

g^ = diag(0,0,0,0,l,-l) . (2.10c) 

With these choices, one sees that 

3 

XI /^P<95 = diag(/xi,-/xi,/X2,-/^2,/^3,-/^3) • (2.11) 
p=i 



Hence, X]p=i A*p Qp assigns eigenvalue +iip to 0p and assigns —fip to the complex con- 
jugates 0*. For the fermions, Y,]>=i Qt ^ diag(/ii, /i2, /^a, ft'^) with 



/ii = |(/ii +/i2 + /^3), h = Uf^i- (^2- (J^s), (2.12a) 
At3 = f(-A*i +A*2 - A^a), /i4 = |(-Aii - A*2 + A*3), (2.12b) 

so the fermion Xj has an effective chemical potential jlj. 
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To derive the correct form of the action to use in a functional integral when chemical 
potentials are present, one could start from the Lagrangian without chemical potentials, 
find the conserved charges of the global symmetry, re-express these in terms of conjugate 
momenta, use this to derive a Hamiltonian path integral representation of the grand 
canonical partition function ( |2.4| ), and then integrate out the conjugate momenta. But 
that is completely unnecessary. It is equivalent, but simpler, to imagine gauging the 
U{lY Abelian global i?-symmetry. The time component of such a fictitious gauge field 
couples to the conserved U (1)^ i?-charge densities, just like the chemical potentials. But 
the time component of a gauge field in a Euclidean functional integral corresponds to i 
times the Minkowski gauge potential. Consequently, in a Euclidean functional integral, 
adding chemical potentials is exactly equivalent to turning on a constant imaginary 
value for the time component of a background gauge field associated with the f/(l)^ 
global i?-symmetry. One ends up with a standard Euclidean functional integral in 
which the Lagrangian contains modified covariant derivatives for the time direction, 

^ D,y - jj,pQp6^o . (2.13) 

It will be convenient for later use to rewrite the Weyl fermions as Majorana 
fermions. Recall that a massless two-component Weyl fermion A, in four dimensions, 
may be converted to a four- component Majorana fermion ijj = (^S)- The corresponding 
terms in the Lagrange density are related via 

A" {Du 5,,o) \^ = li'il?-fx 7075) ^ , (2.14) 

where we have defined ip = (A", Aq.) and 

= (1, ia) , = (1, -ia) , (2.15a) 

lu = (^^ , 75 = 7o7i7273 = Q • (2.15b) 

The Majorana spinors satisfy the condition 

ilj = Cip , (2.16) 

where C = ("p'' ^f^) is the charge conjugation matrix with 612 = —621 = — 1- 

The gauge field may be regarded as a traceless Hermitian matrix. Equivalently, 
it may be expanded as = A^T", with real coefficients {A^} and Hermitian color 
generators T"" which satisfy 

[T\T^]=ir^''T'' , (2.17a) 

and 
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tr {T'^T'') = ^S"'' . (2.17b) 

With this choice, the structure constants /"''^ are real and completely antisymmetric. 

The scalar and fermion fields may similarly be expanded in the basis of color gener- 
ators, (pp = (ppT"" and ipi = ipfT"'. The coefficient is a four-component Grassmann- 
valued spinor. The conjugate spinor is not independent, but is related to ipf via the 
Majorana condition ( |2.16| ). The coefficients 0^ satisfy the reality condition ( ^.6] ) (for 
each a). It will be convenient to introduce explicitly independent real scalar fields via 
(pp = {Xp + iYp)/\^ and to assemble these into a multiplet, 

<^ = {X^,Y,,X2,Y2,X,,Y,) . (2.18) 

We will use a capital Latin index to denote components of this vector, so $a = 
with A = 1, ■ ■ ■ ,6. 

We are finally ready to write down the Af = 4 SYM Lagrange density p5|. 



in Euclidean signature with the addition of chemical potentials and spatial curvature 
induced mass terms: 



£ =tr |i(F^,)2 + {D,Xp - ii2p5,flYpy + {D,Yp + i^XpS^^^Xpf + R-^ ($, 

+ #,(^ - fi, 7075)^^. + \ g^iii^A. ^b]? -gA [{oi%Xp + il3t^i^Y,), I , (2.19) 



where F^^ = d^A^, — d^A^ + ig[A^, A^] and = + ig[A^, ■ ]. The indices run over 
A,B = 1,- ■ -,6, p,q = 1,- ■ -,3, i,j = 1,- ■ ■, 4, and a = 1, ■ ■ -, A^"^ — 1 (with implied sums 
over all indices). The 4x4 matrices and (3'^ satisfy the relations 

{aP, a'^} = -2 6^^ U^, , (3'^} = -2 6^'^ U^, , K, = , (2.20) 

and explicit forms can be given as 

1 / 2 f ^ —<^3\ 3 f i<^2 



= ( .° ?1 , = f ° , = (~T ° 1 ■ (2-21b) 



^1(72 / \~^o / \ ia2y 

The partition function for the grand canonical ensemble has the functional integral 
representation 

VA^ Vtpi V^A ^ , (2.22) 



where the integral is over fields taking values on (of radius R) x (of radius 
/?), with the gauge and scalar fields periodic and the Grassmann- valued fermion fields 
anti-periodic on the Euclidean time circle. 
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2.3 Gauge Invariant Partition Function 



We want to compute the partition function ( p.22|) in the free-field limit of the the- 



ory, while preserving the projection onto gauge invariant states which is embodied in 
the functional integral over Aq. We sketch the calculation here, and leave details to 
Appendix 0. 

Taking a naive g ^ limit of the Lagrangian ( p.l9| ) is the wrong procedure; this 



ignores the special role of Aq as a Lagrange multiplier enforcing Gauss' law and, as 
detailed below, would lead to an ill-defined Gaussian integral due to a zero-mode in 
the covariance operator for Aq. The physically appropriate limit is obtained by sending 
g ^ after separating and rescaling the constant mode of Aq, 

Ao{x) Ao{x) + a/g, (2.23) 

with Ao{x) constrained to have a vanishing spacetime integral, and a an arbitrary 
traceless Nc x Nc Hermitian constant matrix. The net effect is that the g = functional 
integral includes configurations in which parallel transporters around the time circle, 
U = V {e^^ -^0 '^^''^0^^ cover the entire gauge group. It is the resulting integration over 
U which implements the projection onto gauge invariant states. After suitably fixing 
a gauge, the remaining functional integral is a well-defined Gaussian integral over all 
degrees of freedom other than a, the constant mode of Aq, and a non-Gaussian integral 
over the single matrix a. 

The Gaussian integral leads to functional determinants of the covariance (or small 
fiuctuation) operators defined by the quadratic part of the action. For the scalar fields, 
this operator involves the sum of the square of a covariant time derivative, which 
depends on a and the chemical potentials, a spatial Laplacian, and the curvature- 
induced mass term. The gauge field covariance operator and the square of the fermion 
operator have completely analogous forms, but with the relevant Laplacian acting on 
vector or spinor fields, respectively, and without any mass term. Eigenf unctions of 
these operators are exponentials in time, e*"^**, multiplied by scalar, vector, or spinor 
S^-spherical harmonics. The frequencies ujk are quantized Matsubara frequencies, and 
equal 2kTT / (3 for bosons and {2k+l)TT / (3 for fermions. The eigenvalues of the Laplacian 
on the three-sphere may be regarded as squares of discrete spatial momenta, and scale 
as 1/-R^, where R is the radius of the sphere. We will generally set i? = 1, and assume 
all momenta and energies are measured in the units of 1/ R. 

The eigenvalues and associated degeneracies of the spatial part of the small fiuctua- 
tion operator, equal to the spatial Laplacian plus the mass term (for scalars), are shown 
in Table ||. Representations of the SO {A) isometry group are labeled by the integer h 
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field 


eigenvalue 


degeneracy 


transverse vector 


Ax 


{h + lf 


2h {h + 2) 


longitudinal vector 


VF 


h{h + 2) 


{h + lf 


real scalar 


^0 


h{h + 2) 




conformal scalar 


$A 


(h + lf 


{h + iy 


Majorana spinor 






h{h + l) 



Table 2: Eigenvalues (in units of and degeneracies of the spatial part of the small fluctuation 

operator for scalar, spinor, and vector fields. The relevant operator is just the Laplacian on for 
all fields except conformally coupled scalars, where it includes a shift by due to the curvature 

induced mass term. Representations of the 5*0(4) isometry group are labeled by /i = 0, f , 2, • • • , except 
for the longitudinal vector VF, where h starts from 1. 



which runs over all non- negative values; h/R may be regarded as a discrete spatial mo- 
mentum.^ The degeneracy factors are the dimension of the representations at each value 
of h. Spatial components of the gauge fields have been decomposed into a transverse 
{i.e., divergenceless) vector field A_l and a gradient of a scalar, V-F. The longitudinal 
vector V-F vanishes for constant F, and therefore h starts from 1 in this case. The 
spatial small fiuctuation operator for the conformally coupled scalar fields $^ includes 
their curvature induced mass term, so their eigenvalues are h{h^-2) + 1 = (/i+l)^. The 
time component of the gauge field, Aq, is a scalar with respect to the SO (4) isometry 
and has no mass. Note that only Aq has a zero eigenvalue and its associated eigen- 
function is constant on 5*^. The lack of zero modes for divergenceless vector and spinor 
fields is due to the topology of S^, which cannot support constant tensors of rank higher 
than zero. 

The complete small fiuctuation eigenvalues equal these spatial eigenvalues plus the 
square of a Matsubara frequency shifted by an amount proportional to a difference 
of eigenvalues of the matrix a and (for scalars and fermions) by i times a chemical 
potential. Details are given in Appendix ^ For the scalar fields the complete 
small fiuctuation eigenvalues are 

{uJk + qm-qn±il^p)^ + {h + l)\ (2.24) 

where are the real eigenvalues of a, uj^ = 2'Kk/(3, m,n = 1, ■ • • , Nf., and h = 
0, 1, ■ ■ ■ . Note that if the magnitude of the chemical potential /ip exceeds unity (times 



1/R) then, for m = n and h = k = 0, the real part of the eigenvalue ( |2.24|) becomes 



^For the more familiar case of spherical harmonics {i^m} on S'^, the non-negative integer / labels 
the SO{3) irreducible representation whose Laplacian eigenvalue is Z (Z + 1), while the integer m labels 
the basis vectors of the irreducible representation space whose dimension is 2/ + 1. 
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negative — which means that the Gaussian integral fails to be well-defined. This 
can be seen directly from the Lagrangian ( p.l9|) , which shows that non-zero chemical 



potentials act like a negative mass squared for static diagonal components of the scalar 
field ^A- If the chemical potential exceeds the 1/R curvature induced mass, then the 
Euclidean action becomes unbounded below and the partition function ceases to exist. 
So \fip\ = 1/R represents a boundary of the phase diagram of the free theory. In 
the remainder of this section we assume all chemical potentials are less than 1/R (in 
magnitude) . 

The small fiuctuation operator for the static (zero frequency) component of Aq 
has no temporal contribution and is just the scalar Laplacian on S^. The presence 
of a zero eigenvalue of the small fiuctuation operator for Aq illustrates the necessity 
for separating the constant mode of Aq from the other degrees of freedom, as done in 
( |2.23| ). The contributions from the non-constant part of Aq and the longitudinal part 



of the spatial gauge field, VF, end up canceling the contributions from gauge-fixing 
ghosts (which are also scalars with respect to 5*0(4)). The logarithms of the resulting 
functional determinants involve a sum over all Matsubara frequencies and a sum over 
the discrete momentum h labeling 5*^ spherical harmonics. As shown in Appendix 0, 
the Matsubara sum may be performed explicitly and the result for the logarithm of the 
Gaussian integral may be cast in the form of an effective action for the single SU{Nc) 
matrix U = exp(i/?a), 

s,s{u) = -J2- {^^(^") + i-^r^'M^n } Mun ^^u^n - 1] , (2.25) 

n=l 

where x = = e~^l'^ and we have defined the "single particle" partition functions: 

zb{x) = zs{x) + zy{x) , (2.26a) 
, , 6a;2 — 2x3 

zy{x) ^ ^^^--^ , (2.26b) 



X + x^ 



— — - (x^i + x-^^ + x^2 + x-^2 + x^'' + x'^''\ , (2.26c) 
1 — x)'^ 

2x3/2 ^ ^ ^ ^ ^ ^ 

zf{x) = (x2^i + x^aA^i) u-^i^2 ^ 2,-2M2^ u^i^z ^ x~^^^\ . (2.26d) 

(1 — x)'^ 



To obtain the grand canonical partition function, we must integrate over the re- 
maining single matrix a, or equivalently over the group element ?7, 

Z{x) = j dU exp[-S'eff(t/)] . (2.27) 
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The required measure dU is Haar measure on the group SU (Nc). Even though we have 
derived the matrix model specifically for TV = 4 SU{Nc) SYM theory, 



the model dependence is only in the specific field content and group character. The 
generalization to other gauge theories in the zero coupling limit is straightforward. 

2.4 Phase Structure of the Free Theory 

The reduced theory ( p^.27D is a single matrix model. In the large limit, this can 
be solved in a manner similar to the Gross- Witten model [§. The required analysis is 
a straightforward generalization of the zero chemical potential case discussed in Refs. 
[^, Hence we will only sketch the procedure; interested readers should refer to 
Refs. 0, 1 for details.*^ 

After rewriting the integration measure in terms of the eigenvalues of U, and in- 
troducing the eigenvalue distribution function, p{6), one arrives at the following form 
of the effective action for the eigenvalue distribution,^ 

oo 

SeE[p]=N^Y.^-Pn, (2.28) 

n=l 

where pn = J^^d6 p{6) cos{n6) are the Fourier series coefficients of the eigenvalue dis- 
tribution, and^ 

K = ^ {1 - [zBixn + (-1)"+' zp{x-)] } . (2.29) 

The factor of A*"^ multiplying the effective action ( p.28| ) implies that the minimum 
value of this effective action determines the leading large Nc behavior of the free energy 
F = —(3~^ In Z, with fluctuations in the eigenvalue density only generating subleading 



^The analysis in Refs. ||, 0| was carried out with SU{Nc) and U{Nc) gauge groups, respectively. 
But the difference between the gauge groups is negligible in the large Nc limit, and both treatments 



yield the same action (2.28) 



^The eigenvalue distribution p{d) must be non-negative and satisfy the normalization condition 
J^^d9p{9) = 1. The positivity condition implies constraints on the Fourier coefficients {pn} which 
limit how large any particular coefficient can grow. Because of this boundary on the space of allowable 



Pn, the effective action (2.28|) remains bounded below, with a well-defined minimum, even when the 



coefficients Vn are not all positive. This boundary is irrelevant in the disordered phase where all Vn 
are positive and the minimum of Scs lies at p„ — (for n > 0). When one or more of the Vn are 
negative, then the minimum lies on the boundary and its presence is essential. 

^The 1/n term in Vn comes from the Haar measure dU. Expressed in terms of the eigenvalues of 
the matrix U, this generates a Van der Monde determinant which may be written as a contribution of 
-N^ de dO' p{e) p{9') ln|2sin(f-|^)| to S'off[p]. Inserting the identity In 1 2 sin f| = - ^ cosnx, 

and assuming that p{9) is an even function of 9, leads to the form J2^=i n Pn- 
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0{N^) contributions to the free energy. As Nc oo, 



F = min ^ + O(iVO) . (2.30) 
{p} P 



In the effective action ( p.28D each positive Vn acts as a repulsive potential for 



the eigenvalue distribution, while negative Ki's act as attractive potentials. Thus the 
condition 

ZBix"") + < 1 (2.31) 

ensures that all Vn are positive and implies that the system is in the phase where 
eigenvalues repel and the uniform distribution p{9) = ^ characterizes the equilibrium 
state. As long as all | less than one, which is required for the existence of the 

grand canonical ensemble, it is easy to see that our modified single particle partition 
functions (|2.26|) increase monotonically with x. Hence, the n = 1 term in the inequality 



(|2.31|) gives the most stringent condition. Therefore, the curve in the fi-T diagram 
determined by the threshold condition 

zb{x) + zpix) = 1 (2.32) 

gives the boundary of the phase in which the eigenvalue distribution is uniform and 
tr(f/"') has vanishing expectation value (for all n ^ 0). 

In the free-field and large Nc limit, the Polyakov loop expectation value is just pi: 

-l-tTv(^e''^o'^-''Ao^^ = (^^jre'^''^ = (^^JrU^ = J dO p{e) e'^ = pi . (2.33) 

So in the disordered phase, where the eigenvalues of U are uniformly distributed on 
the unit circle, the Polyakov loop expectation vanishes. In the ordered phase, where 
the eigenvalue distribution is non-uniform, the Polyakov loop will have a non-zero 
expectation value. The Polyakov loop transforms non-trivially under the Z^r^ center 
of the gauge group, and its expectation value is an order parameter for the realization 
of this symmetry. Thus the Z^v^ symmetry is unbroken in the disordered phase, and 
spontaneously broken in the ordered phase. As emphasized in the Introduction, it is this 
behavior of the Polyakov loop and the associated realization of the ZAr^-symmetry which 
motivates calling the disordered phase "confining" and the ordered phase "deconfined" 



^The expectation value of the Polyakov loop may be interpreted as exp[— /3Ai^] where AF is the 
free energy difference between equilibrium states in which one has, or has not, added a fundamental 
representation static test quark. In the confining phase this free energy difference is infinite and 
the expectation value of the Polyakov loop vanishes, while in the deconfined phase the free energy 
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Figure 2: The confincment/deconfinement phase transition hne in the fi-T plane for the free theory. 
Both yLt and T are measured in the units of 1/ R. The values of (/xi, /i2, /is) for each line are as follows. 
Solid line: (/i, 0, 0), dotted: (/i, /i/2, /i/2), dash-dotted: (/i, /i, 0), dashed: /i). 



Returning to equation ( |2.32D , it is straightforward to plot where, in the /i-T plane, 
solutions to this condition lie. Fig. |^ shows the resulting curves in four representative 
cases where (/xi, /i2, /^s) = (At,0,0), (/x, /i, 0), (/x, or (/i, /i/2, /i/2). Since there are 
three independent chemical potentials, the full phase diagram is four dimensional, but 
the slices shown in Fig. |] illustrate the general behavior. When all chemical potentials 
are zero, the phase transition line reaches the T-axis at T 0.38/i? agreeing, as it 
should, with Ref. 0. If |/i| denotes the maximum magnitude of the three chemical 
potentials then, as discussed earlier, |/i| = 1/R is a boundary of the free field phase 
diagram. In the vicinity of zero temperature, the relevant terms in Eq. (|2.32| ) are 

g-/3(l-Mi) _^ g-/3(l-M2) _|_ g-/3(l-At3) _|_ 2 g-5/3(3-Mi-/t2-/t3) _ ^ (2.34) 

where the first three exponentials are from zs{x), the last one is from zf{x), and 
the radius of the sphere R is set to unity. From this expression, one sees that the 

difference is finite and the expectation value is non-zero. Hence the Polyakov loop may be viewed as 
a confinement/deconfinement order parameter. But strictly speaking, the Polyakov loop expectation 
value is defined via cluster decomposition from the large distance limit of the Polyakov loop two-point 
function. It is the two-point function which, in infinite volume, embodies the operational definition of 
confinement in terms of the free energy needed to separate a test quark and antiquark to infinity. 
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maximal chemical potential must approach 1 as T ^ 0. Thus, the phase transition line 
necessarily ends at the boundary = 1 when the temperature falls to zero. The slope 
of the transition line at zero temperature depends on how many chemical potentials 
approach unity. It is easy to see that the limiting slope of the transition line is zero if 
a single chemical potential is turned on. The limiting slope is — In 2 when two equal 
chemical potentials are turned on, and — In 4 with three equal potentials. 

2.5 Order of the Phase Transition 

In the confining phase, the density of eigenvalues of the matrix U is constant and all 
(non-trivial) Fourier coefficients of p{6) vanish. The minimum value of the effective 
action ( p^.28| ) vanishes, and hence the free energy in this phase is 0{N^), not 0{N^), 
as A'c — ^ oo. In the deconfined phase there is a non-constant density of eigenvalues, 
the minimum value of the effective action (|2.28|) is non-zero, and the free energy is 



0{N^). In other words, limAr^^oo F/N^ is non-zero in the deconfined phase but vanishes 
identically in the confined phase. The free energy must be continuous across any 
phase transition (but its derivatives need not be), so the coefficient of the 0{N^) part 
of the free energy must vanish as one approaches the phase transition line from the 
deconfined side. If it vanishes linearly (with temperature or chemical potential) then 
first derivatives of the free energy will be discontinuous and the Nc = oo phase transition 
is first order. If the 0{N^) free energy vanishes faster than linearly as the transition 
line is approached, then the = oo phase transition is continuous. 

The phase transition line is determined by the condition zb{x) + zf{x) = 1, and 
the left-hand side is a monotonically increasing function of temperature. Therefore 
as we approach from the deconfined side, we can analyze the local behavior near the 
phase transition line by expanding the solution to the matrix model in powers of = 
zb{x) + zf{x) — 1, with e real. Such an expansion was carried out in Ref. ^ for the 
case of zero chemical potentials, with the result that (for > 0) 

lim ^ = -T + (^(a- (2.35) 

(Additional higher order terms are also obtained in Ref. 0.) The analysis is indepen- 
dent of whether the single particle partition functions contain chemical potentials, so 
the above result is equally valid in our case with chemical potentials. 

To see if the phase transition is first order, one must merely determine how in 
the result (|2.35|) depends on T — Tc or /i — /ic. The function = zb{x) + zp{x) — 1 
depends analytically on T and /i and there is no reason for its derivatives with respect 
to T or /X to vanish on the line where crosses zero. It is straightforward to check 
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numerically that this is, in fact, the case. One finds that vanishes linearly with 
T — Tc and /i — /Xc as one approaches a point (Tc, /Xc) on the phase transition line from 
the deconfined phase. Hence the rescaled large N^. free energy, limAr^^oo -^/^c 5 has a 
discontinuous first derivative as one crosses the phase transition line, showing that the 
large N^. confinement /deconfinement phase transition, at zero coupling, is first order. 

First order phase transitions are normally robust phenomena with regard to per- 
turbations, such as a change in the coupling in the underlying theory. So one might 
expect that a first order confinement /deconfinement transition in the zero coupling 
limit of the theory would imply that the transition must remain first order for some 
non-zero range of couplings. In the case at hand, however, the situation is more subtle. 

At a first order transition, there are multiple equilibrium states so that, as the 
transition is approached, the limit depends on the direction of the approach. At a typi- 
cal first order transition, the coexisting equilibrium states are separated by free energy 
barriers, and it is the presence of these free energy barriers which lead to characteristic 
phenomena associated with first order transitions such as superheating or supercooling. 
But m. Nc = 00 gauge theories at zero coupling, such as the J\f = 4 SYM theory under 
discussion, the free energy at phase coexistence does not have isolated minima, but 
rather has a fiat direction. This may be seen directly in the effective action (p. 28 ) for 



the density of eigenvalues. The coefficient Vi vanishes at the deconfinement transition, 
so minimizing the effective action leaves the Fourier coefficient pi completely undeter- 
mined. Approaching the phase transition line from the confined or deconfined phases 
leads, at the transition, to the equilibrium states with minimal or maximal values of 
pi, respectively. But these states are not separated by any free energy barrier. 

The lack of a free energy barrier separating the coexisting equilibrium states (in 
the zero coupling theory) means that a small perturbation, such as turning on an 
infinitesimal non-zero coupling, can have a large effect. A non-zero coupling can lift 
the flat direction and produce either a flrst order transition (with a barrier separating 
coexisting states) or a second order transition, depending on the sign of the (pi)^ term 
which is induced in the effective action for the density of eigenvalues. As discussed in 



^'^For example, with only one chemical potential turned on, (TcMc) = (0.35,0.53) (measured in 
the units of 1/R) is a point on the phase transition line and ~ — 2.4(T — Tc) — 0.31 (/i — /ic). 
When all three chemical potentials are equal, (Tc, /ic) — (0.35, 0.32) lies on the transition line and 
^ ^ -2.5 (T-Tc)- 0.48 i^i-^ic). 

^^Because of this, some might quarrel with calling this a first order transition, even through the free 
energy shows a kink (in temperature or chemical potential) as one crosses the transition line. Following 
the analysis of Ref. 0, one may show that if the transition line is approached from the deconfined 
phase, then the Polyakov loop expectation value pi has a limiting value of 1/2, independent of the 
chemical potential. 
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greater detail in Ref . , a three- loop calculation is needed to determine this sign. The 
required three-loop calculation has not yet been done (even for zero chemical potential), 
so the order of the confinement /deconfinement transition in A/" = 4 super Yang- Mills 
theory at small but non-zero coupling is currently unknown. 

2.6 Canonical vs. Grand Canonical 



To understand how our results compare with the picture of Ref. |^ , one must convert 
from the grand canonical to the canonical ensemble (or vice- versa). Starting from the 
grand canonical free energy F(/5, /ip) = — In Z(/?, /Xp), the -R-charges are given by 
Qp = —■§^, and the thermodynamic potential of the canonical ensemble is A{P, Qp) = 
^ + '^pl^pQp (with the chemical potentials now viewed as functions of the charges). 
Given the canonical ensemble potential as function of temperature and charge, the 
chemical potentials are given by /ip = . 

Within the deconfined phase, the free energy and the i?-charges are both of order 
iV^, so the rescaled charges 



are appropriate observables in the large Nc limit. Ref. ||2^ considered thermodynamics 
with a fixed non-zero value for the total rescaled charge q = 'YlipQv 

Consider, for simplicity, the case of equal chemical potentials. In the interior of 
the deconfined phase, there is a one-to-one mapping between the chemical potential /i 
and the rescaled charge q. But within the confined phase, the i?-charges (as well as 
the free energy) are 0(1), so the rescaled charge q vanishes identically. At any fixed 
temperature T below the = confinement transition temperature, if one increases 
the chemical potential starting from zero then q remains identically zero until one 
reaches the confinement /deconfinement transition at /i = fJ^dT)- Since the transition 
(in the free theory) is first order, the rescaled charge q jumps discontinuously across 
the transition to some non-zero value qc{T), which is the minimal value of q (at the 
given temperature) within the deconfined phase. As the chemical potential is further 
increased, the charge continues to grow, eventually diverging at the edge of the phase 
diagram {i.e., as /i — >• 1). 

From this description, it might seem that it is impossible to obtain an equilibrium 
state with < g < qc{T) — but this is wrong. The essential point is that first order 
transitions always involve phase coexistence. The q = qdT) and g = statistical 



^^For pure Yang-Mills theory on a sphere, the corresponding three-loop calculation has been per- 
formed 1^ , and in this case the large Nc confinement transition remains first order at small but non-zero 
coupling. 
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Figure 3: Schematic phase diagram for the free theory in the q-T plane. Within the phase coexistence 
region, the chemical potential is independent of the charge q, and equal to its value iJ-c{T) on the 
confinement/deconfinement transition line. Inside the deconfined phase, the chemical potential (at 
fixed T) increases with increasing q. Only the T < Tc segment of the q = axis corresponds to the 
pure deconfined phase. The maximal charge q^^^ = 1/4, and the maximum transition temperature 
ymax ^ 0.38/i?. 



ensembles produced by taking fi — > fic{T) from above and below, respectively, in the 
grand canonical ensemble are extremal equilibrium states. But any statistical mixture 
of these two states is also a valid equilibrium state in the Nc —>■ oo limit. 

Therefore, the canonical description of the phase diagram shown in Fig. |^ will have a 
deconfined phase for q > qdT) in which the susceptibility ^ = ^ is positive (so the 
chemical potential increases monotonically with q), together with a phase coexistence 
region for < g < qdT) in which the the chemical potential n = limTv^^oo ]^ is 
independent of q. Within the phase coexistence region the (rescaled, large Nc limit 
of the) thermodynamic potential, limN^^^oo A / , is simply equal to fi(T)q. Only at 
g = will one have the pure confined phase. This is illustrated in Fig. § 

Although the authors of Ref. did not explicitly evalute the chemical potential 



produced by their thermpdynamic potential, or notice its independence on q for their 
"small g" extremum, the results of Ref. |^ are completely consistent with the above 
description of phase coexistence in the q-T phase diagram. 



""^•^ Within the truncation used in Ref. , the thermodynamic potential in the confined phase is given 
by f3A{f3, q) = miup [A^^ — Scs[p)\ , where Sch{p) is their Eq. (3.19). Evaluating this gives a result 
proportional to q, and hence a chemical potential independent of q. These authors also constructed 
a phcnomenological model based on a deformation of the effective action which can mimic the effect 
of turning on a weak gauge coupling. For this generalization, assuming the transition remains first 
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3. Weak Coupling, High Temperature Effective Theory 



We now consider the theory in the high temperature regime, T ^ with a small 
non-zero coupling. The relevant expansion parameter will be l/(Ti?)^, so the high 
temperature regime may also be viewed as the regime of large spatial radius at a fixed 
temperature. Consequently, for observables which are primarily sensitive to the scale 
T, the curvature of the will produce only small corrections to fiat space results. 

Let A = g^Nc denote the usual 't Hooft coupling, which is held fixed as Nc ^ oo. 
Interesting phenomena in the phase diagram will be found to occur when l/(Ti?)^ and 
/i^/T^ are both of order A, so we will focus on this region of parameter space. 

At non-zero temperature, one may always decompose fields into Fourier series run- 
ning over all Matsubara frequencies, and thereby view every four-dimensional field as 
an infinite tower of three dimensional fields. At high temperature, where the 1/T cir- 
cumference of the periodic time circle is small compared to all other scales, the non-zero 
Matsubara frequency modes behave (with respect to physics on length scales large com- 
pared to 1/T) like fields describing very heavy excitations with 0(T) mass. As T — oo, 
the non-static heavy modes decouple and one may integrate them out leaving only 
the light zero frequency modes. In this way, one obtains a three dimensional effective 



theory of the static modes |]28[. Convenient techniques for systematically constructing 
this effective theory, using dimensional continuation to regulate both long and short 
distance behavior, were described by Braaten and Nieto [EI 



There are three scales which are relevant for equilibrium thermodynamics in weakly 
coupled high temperature non-Abelian gauge theories at zero chemical potential: the 
temperature T, the "electric mass" (or inverse Debye screening length) which is C(v^T), 
and the "magnetic mass" (or inverse correlation length of the static gauge field) which 
is 0{XT). One may construct a dimensionally reduced "electric" effective theory, valid 
on length scales large compared to 1/T, by integrating out non-static fluctuations. One 
may then construct a further reduced "magnetic" effective theory, valid on length scales 
large compared to (VXT)^^, by also integrating out static fluctuations on the Debye 
screening scale. For QCD, these two effective theories are commonly referred to as 
EQCD3 and MQCD3, respectively 13 . 

For our case of A/" = 4 super- Yang-Mills theory with non-zero chemical potentials, 
it is the "electric" effective theory which will be of interest. This effective theory, 
which we will call "ESYM3", will contain a three dimensional gauge fleld A, the static 



order, one again finds that A is linear and the chemical potential is independent of q in the region 
of the phase diagram where the thermodynamic potential is minimized by the stable small q solution 
(denoted solution / in (22|). Hence this region again corresponds to phase coexistence, and the phase 
diagram resulting from this model resembles Fig. above. 
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component of Aq, which will appear as an adjoint scalar field, and the static components 
of the original scalars The fermions, having no zero- frequency components due to 
their antiperiodicity, will be completely integrated out. Operationally, one builds the 
effective theory by considering all local operators which may be constructed from these 
fields and are consistent with the symmetries of the theory. Up to any given order in 
perturbation theory, only a finite number of operators with sufficiently low dimensions 
are needed. The required coefficients are determined by matching the results for a 
minimal set of observables which may be evaluated in both the effective and underlying 
theories [^, [S^ . 



The operators in ESYM3 include the identity operator, the usual gauge invariant 
derivative terms, quadratic mass terms for scalar fields, and scalar interaction terms. 
It is convenient to use rescaled fields in the effective theory. 



A = ^/Z,TA, Ao=^/Z2TAo, ^A = VZ3T^A, (3.1) 

with Zi, Z2, and Z^ dimensionless wavefunction renormalization factors, and and 
now having canonical dimension 1/2. The wavefunction renormalization factors 
may be chosen so that the Lagrangian density of the effective theory has the form 

/^ESYMa = / + tr + (Aio)' + MlAl+ (A$a)' + ml + V{Ao, ^a) , 

(3.2) 

where Di = di + igslAi, ■] and Fij = [Di, Dj], with g-^ = gVT the dimensionful gauge 
coupling appropriate for a three-dimensional theory. 

The coefficient of the identity operator, /, is the 3-d effective theory version of a 
cosmo logical constant, and represents the contribution to the free energy density from 
modes which have been integrated out. Because we are working at high temperature 
(or large volume) and only integrating out heavy modes with short 0{1/T) correla- 
tion lengths, the computation of the free energy density due to heavy modes may be 
performed entirely in flat space, treating the curvature and chemical potential induced 
scalar mass terms as small perturbations. Corrections to this approximation will vanish 
exponentially fast in TR (or faster than any power of A, since we are assuming that 
TR is of order l/y/X). Details of the calculation are described in appendix p.l| . The 
result, through order A, is 

/._^|2.^-3a-^+x:|f+i:|I+o(a4' p-^) 

p=l i=l 

The mass parameters m\ of the scalar flelds $a receive tree-level contributions from 
the curvature coupling and the chemical potentials, plus 0{\/XT) one-loop thermal 
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corrections, so that 

m\ = R-^ - ^i\ + 5m'^{T), (3.4) 
where we introduce (for later convenience) 

/i2 , A = 3 or 4; (3.5) 



AH , 


A 


= 1 


or 


2 




A 


= 3 


or 


4 


/^3, 


A 


= 5 


or 


6 



Because we will be interested in the regime where the tree-level contributions nearly 
cancel, including the thermal mass correction 5m? {T) is essential. This contribution, 
as well as the one-loop scalar wavefunction renormalization (which will also be needed) 
may be obtained by matching the two-point scalar correlator in the original theory 
at zero frequency and low spatial momentum with the corresponding correlator in the 
effective theory. The computation is shown in Appendix [B.2| . The result for the thermal 



mass (also derived in Refs. |3T|, |3^, Q) is 



W(T) [A + C(A')] . (3.6) 

One-loop fluctuations on the scale T also produce an 0{-\/~XT) mass for Aq. The 
resulting Debye mass parameter is [pT| 



Ml = T'^[2\ + 0{\^)\. (3.7) 

The interaction potential V{Aq^ ^ a) contains both tree-level and fluctuation in- 
duced thermal contributions, 

V^(io, ^a) = Krcc(io, ~^a) + 5V{A,, ^a) ■ (3.8) 

The tree-level part equals the scalar interaction terms of the original four- dimensional 
action, re-expressed in terms of the rescaled fields (|3.1|) (with the wavefunction renor- 
malization factors equal to one at lowest order), 

Krce(io, ^a) = tr {2(73 /ip ([^,^p] Yp) + qI {i[A^, ^a]? + \gl (^[$A, ^B\f] . (3.9) 

It is straightforward to obtain one- loop contributions to 5V that are quartic in The 
coefficient of the quartic term can be computed by evaluating the one-loop heavy-mode 
contribution to the four-point correlator of the scalar fields with vanishing external 
momenta. Once again, because we are interested in effects due to heavy modes with 
correlation lengths small compared to i?, these contributions may be evaluated in fiat 



- 25 - 



space, ignoring both curvature and chemical potential corrections. The required calcu- 
lation is performed in appendix p.3| , with the result 



5Vr 



quartic 



(3.10) 



where the trace is in the adjoint representation, with {^^^)ac = i^^/""^. Note that no 
terms cubic in $^ can be generated, since they would not respect the U (1)^ subgroup of 
the SU{A) i?-symmetry, which is preserved even in the presence of chemical potentials. 

One may, in principle, continue in a similar manner and obtain the coefficients 
of arbitrarily higher dimensional operators in the effective action. However, in the 
special case where the scalar fields take values in flat directions (which minimize the 
tree potential), one may employ a more elegant background field method. This method 
not only confirms the coefficients of the results ( |3.(iD and ( ^.lUD , it also yields all higher 
order terms, in powers of $, in the one-loop contribution to 6V{^a)- This is discussed 
in Appendix and the result is 



ivr^ tr 



(=3 



+ J]8(l-4-'+2)^^yi^C(2/-3) 



(2/)!! 



„2 „ _ \ « 



7r2T2 



(3.11) 



These adjoint representation traces may be written more explicitly in terms of the 
eigenvalues {A^} of (viewed as an Nc x Nc matrix). It will prove convenient to 
introduce dimensionless rescaled eigenvalues. 



93 \m 
A 1 



(3.12) 



in terms of which 



tr 



^2 

^2-p2 Z_^^A ^A 



(3.13) 



Note that when the rescaled eigenvalues {p™} (and their differences) are parametrically 
of order one, then all terms in the series ( p.ll| ) are comparable in size. The potential 
(|3.11|) will play an essential role in the next section, where its behavior will be examined 
in more detail. 
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4. High Temperature Thermodynamics 



4.1 Deconfined Plasma Phase 



The effective scalar masses m\ [given by Eqs. (^) and ( |3.(i| )] are positive as long as 



the chemical potentials are sufficiently small, 



fil<R~^ + XT\ 



(4.1) 



Assuming this is the case, the trivial configuration of vanishing scalar fields. 



— 







(4.2) 



is a local minimum of the effective potential. The conditions under which this is the 
global minimum will be examined in the next subsection. Because of the non-zero Debye 
mass M|, [given by Eq. (|3.7|) ], the static mode Aq has perturbatively small fluctuations 
in the weak coupling, high temperature regime. Hence, the Polyakov loop expectation 
value (trf/(a;)) = (|tr e*^^"^"^"^-*) is non-zero. So the trivial configuration (|4.2| ) describes 
a deconfined plasma phase in which the center symmetry is spontaneously broken while 
the U{1)^ i?-symmetry (left invariant by the chemical potentials) is unbroken. 

The free energy F = —T In Zesyms , where ^esyms is the partition function of the 
three-dimensional effective theory ( p.2| ). In the deconfined plasma phase the order A° 
and contributions to /3F arise solely from the heavy modes which were integrated out 
to produce the effective theory. Hence, these contributions are completely contained 
in the coefficient of the unit operator, /, given in Eq. ( |3.3| ). The next contribution is 
of order A^''^ (not A^), and comes from static fluctuations of the scalar fields and 
Aq. In the full theory, evaluation of the order A^^^ contribution requires an infinite 
resummation of ring diagrams. But this becomes an easy one-loop calculation using 
the high temperature effective theory. 

Because the spatial gauge fields Ai are massless, their one-loop contribution van- 
ishes (using dimensional continuation). The effective scalar fields Aq and $a each make 
a one-loop contribution to the free energy density of the form 



where the integral is evaluated using dimensional continuation, and m is the appropriate 
mass parameter for the field. So the entire (9(A^/^) contribution to the free energy 




(4.3) 



i-^Explicitly, {trU) = Nc[l + i'^^w^)'^^ + C(A^)] in the deconfined plasma phase. This follows 
from the pure Yang-Mills result |Q after accounting for the differing Debye mass of TV = 4 SYM. 
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density is 



127r 



D 

[mi + 5^ m^) . (4.4) 



A=l 



Inserting the explicit expressions (|3.4|) - (|3.7|) for these masses, and adding (T times) 



the coefficient / of the unit operator ( |3.3|) , we obtain the free energy in the deconfined 
plasma phase, through order A^/^, 

r2 r Q 3 ^,,2 4 ^ 

2^2 



-f^ plasma - 12 [ T2i?2 2^2 ^ 



(4.5a) 



with V = 27r^i?3 the spatial volume. Note that this expression is only valid (and real) 
when the chemical potentials satisfy the inequality ( [4.1|) . 



4.2 Is There a Higgs Phase? 

Now consider the situation when the maximal chemical potential exceeds AT^ + i?"^, so 
that some of the effective scalar masses m\ = AT^ + — ^\ become negative. This 
makes the quadratic terms in the ESYM3 scalar potential unstable, so the minimum 
of the effective theory scalar potential can no longer lie at the origin. However, one 
may verify that the sum of the tree ( ^.9] ) and one-loop ( p.lUD quartic contributions to 



the scalar potential is positive definite for all non-zero showing that the quartic 
thermal corrections to the effective potential lift the flat directions which are present 
in the tree-level potential. 

This suggests that thermal corrections to the scalar interactions in the effective 
theory may stabilize the theory when chemical potentials are large, generating a non- 
trivial minimum of the effective scalar potential. Such a non-trivial minimum would 
correspond to a Higgs phase with spontaneously broken i?-symmetry. 

Assume (for the moment) that terms in the scalar potential involving higher than 
quartic powers of the field are unimportant. One can explicitly minimize the quadratic 
plus quartic contributions to the ESYM3 scalar potential. At the minimum one is 
balancing an unstable tree-level quadratic term against a one-loop quartic term, so the 
rescaled eigenvalues ( |3.12| ) of the scalar fields $a are of order 



{PA? - ^ • (4.6) 
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The corresponding value of the (truncated) potential, at its minimum, scales as 



rp3 



max 



(- 



m 



(4.7) 



A V AT2 

However, minimizing the truncated scalar potential and adding its value to the coeffi- 



cient of the identity operator (|3.3|) does not produce a valid approximation for the free 
energy and, as we will see, this putative Higgs phase simply does not exist. 

There are two problems with the above scenario suggesting the presence of a 
Higgs phase. First, we neglected higher than quartic terms in the ESYM3 potential 
(|3.11|) . This is only acceptable if (p™)^ <^ 1 which, given the above estimate, requires 
\mW <^ AT^. The curvature and chemical potential contributions to m\ are, by previ- 
ous assumption, individually of order AT^. So neglecting terms higher than quartic in 
the potential requires that m\ be parametrically smaller than its individual pieces. In 
other words, one must restrict attention to a "fine-tuned" region in the phase diagram 
sufficiently close to the surface where an effective scalar mass ffist turns negative. 

Second, and more importantly, the above scenario was based on a classical analysis 



of the action (|3l^ ) of the effective theory — so it neglected the effects of fluctuations 
in the static ESYM3 fields. To see if this matters, let 

$A = {^a) + 6^A , (4.8) 

with the mean value {^a) lying along a flat direction of the tree- level potential. That 
implies that the mean values are simultaneously diagonal, (($A)mn) = {t^T / g^) p'^Smn 
(up to an irrelevant gauge transformation). Inserting the decomposition ( |4.8| ) into 
the tree-level potential ( |3.9| ) leads to non-zero mass corrections for the off-diagonal 
components of the scalar field and Aq fluctuations,^^ 

6VtU^A) = Ml^[\{Ao)mn\' + J2\{S^B)mn\'} + 0(5<l=^) + 0{6^AI) , (4.9) 

m,n B 

with 

Ml^^i.'T'Y.^p^-pAf- (4.10) 

A 

And inserting ( |4.8|) into the scalar kinetic term produces mass terms for the off-diagonal 
components of the static gauge field, tr(A$yi)^ = ^m,n ^mn \ {A)mn\ H • 



^^We have omitted ofF-diagonal terms contained in SVtrcci^ a) which mix S^a with Aq, or mix S^a 
with (5$B (for A ^ B). These terms may be canceled by using the three dimensional version of the 



gauge fixing term (C.2), with gauge fixing parameter ^ = 1. 
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If the (rescaled) diagonal components {p^} are 0{1), then the mass terms (f4.10|) 
for off-diagonal fields induced by the mean values of the diagonal components will be 
of order T^, or large compared to the (9(AT^) thermal contributions due to the non- 
zero frequency modes, as well as the curvature and chemical potential induced mass 
terms [which we have assumed to also be (9(AT^)]. It will be sufficient to focus on this 
regime; if a Higgs phase does exist, then the typical size of will need to be of this 
order. Consequently, if the scalar fields $a have non-trivial mean values lying along 
flat directions of the tree-level potential, then off-diagonal components of the static 
fields will act like heavy degrees of freedom (just like the non-static Matsubara modes) 
and need to be integrated out before considering the dynamics of the remaining "light" 
diagonal modes. 

Doing so is straightforward; the basic ingredient is just the three dimensional loop 
integral ( [4.3|) . The result is a "Higgs-branch" three-dimensional effective theory which 
we will term "HSYM3" and whose Lagrange density has the form 



HSYM3 



tr 



(4.11) 



where the fields are now diagonal, with gauge group [U{l)]^--\ The scalar potential 
V{^a) for the remaining diagonal components is 



V{^a) = f + tr(ml $1) + 5Vflat($A) + 5Kff-diag($. 



(4.12) 



where / is the additive constant (|3.3|) due to the non-zero frequency modes, 5Vflat($A) 
is the one-loop potential ( p. Ill) induced by non-zero frequency fluctuations, and^^ 



ofr-diag('^'A 



-fvr^ tr 



A 



(4.13) 



Writing this explicitly in terms of the rescaled eigenvalues {p™} of ^a, the potential is 



V{p 



m,n 



(4.14) 



^^For configurations in which the scalar field has degenerate eigenvalues (which may naturally 
be grouped together into blocks of equal eigenvalues), it is the block off-diagonal components of the 
static fields which acquire mass, and the block-diagonal components which remain light. For these 
exceptional configurations, the residual gauge group is larger than [(/{l)]^"^^ . 

^^This is just the result (4.3), adapted to the differing masses ( 4.10| ), and multiplied by 8 to account 
for the 8N^ (real) bosonic degrees of freedom — 6 from the scalars and two from the transverse 
components of the gauge field. (The gauge-fixing ghosts cancel the Aq and longitudinal contributions.) 
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where we have defined an rms measure of eigenvalue differences, 

p--^[Y^{p--piy\"'\ (4.15) 

A 

introduced the rescaled mass shift 

^ ^ - 1 = , (4.16) 

and defined the one-dimensional function 

h{v) = -\ + v^-iy + {\r.2)v' + j2^{l- 4-'+^) C(2/-3) {-v^)^ . (4.17) 

The additive constant — ^ comes from /. We have dropped subleading 0(AT^) con- 
tributions to /, as these are parametrically smaller than the 0{T^) terms retained in 
V. 

Inserting the defining series for the zeta function and interchanging summations 
allows one to re-express h{v) in an alternative form which is more convenient for ex- 
amining its global behavior, namely 

oo 

h{v) = -l + v'-lv' + IJ2 i-^y-' [if+vY' - f - Iv'j] . (4.18) 

i=i 

The function h{v) is graphed in Figure |[ For large arguments, h{v) approaches zero 
exponentially rapidly, 

h(v) -^(2vf/^e-^\ (4.19) 

The approach to zero as f — > cxd may be understood as a consequence of supersymmetry. 
In the potential ([4.14|) , the argument p"*" of the function h is the ratio of the mass scale 



Mmn for the off-diagonal fields to (vr times) the temperature. A large ratio is equivalent 
to sending T — while holding this mass scale fixed. (Since we are assuming 1/{TRY ~ 
p^ jT'^ ~ A, this also implies sending _R — > oo and /i 0.) At zero temperature, 
supersymmetry requires that all loop contributions to the effective potential cancel. 
So the cancellation of contributions to the one-loop potential as eigenvalue differences 
become large follows from the cancellation of one-loop contributions as T — > 0. 



^^To derive this asymptotic form, it is convenient to rewrite the infinite sum in (4.18) as a contour 
integral, leading to h[z) = -i + w^ + f lime^o |f c.-&c.{t^z)\{z'^ ^v^f^"^ ~{z^ ^t^fl'^ -^v^i^z^ + 
where the contour C encircles the real axis counterclockwise and the branch cuts of the integrand run 
outward from the branch points to ±zoo. Deforming the contour so that it wraps around the branch 
cuts, producing an integral of the discontinuity across the cuts, leads to the stated asymptotic form. 
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Figure 5: Plots of the scalar potential V, divided by the overall factor c = ^n^N^T^, for the 
eigenvalue configuration p™ — v S\ (—1)™, as a function of v. The four plots correspond to SniA^ — 
0.2, 0.0, —0.05 and —1, respectively. Note that half of the eigenvalues differences p™" are zero 
in this configuration, and each vanishing difference contributes h{0) = to V / {^n'^T^). This 
accounts for the non-zero asymptotic value in case (b). If Srh'j^ is positive (or zero), then the potential 
is monotonically increasing with a unique minimum at the origin. If —1 < drrij^ < 0, then the 
potential is unbounded below, with a local minimum at the origin. If Sm^ < ^1, then the potential 
is monotonically decreasing with a maximum at the origin. 
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As indicated in Eq. ( [4.14|) , the complete potential V{p^) is a sum of values of h{v) 
plus the sum of squares of the eigenvalue differences weighted by the mass shift Srh\. 
Plots of the resulting potential, for a simple eigenvalue distribution and several choices 
of mass shifts, are shown in Figure ^.^^ Since the function h{v) is bounded (for positive 
real v), it is evident that the potential V{p^) will be unbounded below if any mass shift 
6m\ is negative. Since Srh^ oc — fi\, this means that no truly stable equilibrium 
phase can exist if any chemical potential exceeds 1/R. Moreover, the potential V^p"^) 
has no local minimum for any non-zero set of eigenvalues {p^}. Hence there is no Higgs 
phase (not even a metastable Higgs phase) in this weakly coupled, high temperature 
theory. 



4.3 High Density Metastability 

When the largest chemical potential lies in the interval 

R-^ < XT^ + R-^ , (4.20) 

the deconfined plasma phase (corresponding to the origin of the field space) is locally 
stable, but is not the global minimum of the free energy — as there is no global 
minimum. For this range of chemical potential and finite Nc, as we will see, the 
deconfined plasma is metastable with a lifetime which grows exponentially as —>■ oo. 
The onset of instability, at a (maximal) chemical potential of l/i?, is independent of A 
and thus exactly where the free theory becomes ill-defined. This may be understood as 
a consequence of the existence of towers of BPS operators with vanishing anomalous 
dimensions (independent of A) which map to towers of BPS states on with energies 
precisely equal to their i?-charge. Consequently, the Boltzmann sum representation 
of the partition function ceases to converge if any chemical potential exceeds 1/R. 
However, this does not prevent the existence of an arbitrarily long-lived metastable 
phase above this threshold. 

To estimate the decay rate when p > 1/R, consider the behavior of the effective 
potential V{^a) for configurations of the form 

p7 = ^SA,i{Sm,i~l/N,), (4.21) 

for some dimensionless real number v. All but one of the eigenvalues remain near the 
origin in this configuration as the single eigenvalue p\ = v (1—^) is varied. (This 
configuration satisfies the SU{Nc) tracelessness constraint, J2mPA = 0-) Assume that 

^^Fig. 1^ plots the potential when the eigenvalues of one scalar field are non-zero and all have the 
same magnitude, with half of them positive and half negative (so that the field is traceless). 
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Hi is the largest chemical potential. For this class of configurations, the potential V{^a) 
becomes 

V{v) = -l-K^T^Nl + -K^T^N, [ h{v) + \ + 5m\ v^] + 0{N^,) . (4.22) 

The interval ( |4.20| ) corresponds to 6rhf = — /i^)/(AT^) lying between and —1. A 
plot of V{v) resembles case (c) of Fig. ^, with a potential barrier separating the region 
near v = from the bottomless region at large values of v. However, for the eigenvalue 
configuration ( ^.21| ) the height of this potential barrier grows only linearly with N^, 

AV = maxV{v) - F(0) = 0{N, T^) . (4.23) 



As seen from Eq. ( [4.14|) , this potential difference receives contributions from every pair 



of unequal eigenvalues of $i and, for this configuration, only 2(A^c~l) terms involving 
the difference between p\ and the other Nc—1 eigenvalues of $i contribute. 

The deconfined plasma phase corresponding to f = can decay via a thermal 
fiuctuation in which a single eigenvalue makes a large excursion from the origin to 
top of the barrier, uniformly in space. The probability for such a fiuctuation has a 
Boltzmann suppression factor of 

^-VAV ^^-OiN.iTRf)^ (4.24) 

where V = 2it^R^ is the spatial volume. (The usual factor of /3 = 1/T appearing in 
a Boltzmann factor was absorbed in our definition of the scalar potential of the three 
dimensional effective theory.) 

Alternatively, a thermal fiuctuation could nucleate a critical bubble in which a 
single eigenvalue makes a large excursion to values across the barrier with lower free 
energy density, over a sufficiently large spatial region so that bubble subsequently grows. 
The process is characterized by a Euclidean bounce solution |]35| , The required size 



of the critical bubble scales as (-\/AT)~^. The action (in the effective theory) for such a 
solution is 0{Nc/ X^^"^), so the rate for this process will have a Boltzmann suppression 
factor of^° 

g-Sbouncc = ^-©(AfcA-a/a) _ (4.25) 

When 1/{TR) = C(a/A), which we have assumed, then both processes will have 
rates which are exponentially suppressed by an amount which scales as Nc/X^^'^. Which 



process dominates depends on the pure numerical coefficients in the exponents of (f4.24|) 



^°To see this, it is convenient to rescale spatial coordinates x — + y/\/\T, which corresponds to 
measuring distance in units of the Debye length. Then, for a single eigenvalue excursion, the scalar 



part of the action ( 4.11 ) reduces to a dimensionless action of the form J (Py [^(Vp)^ + h{p)] times an 



overall factor of Nc/X'^^^- Therefore, the action of the bounce solution must scale in the same way. 
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and ([4.25|) [which we have not evaluated] together with the value of \/\TR?^ Conse- 
quently, despite the fact that no true equilibrium state exists when /i > 1/-R, as long 
as the maximal chemical potential is within the interval ([4.20|) , the deconfined plasma 
phase will be metastable with a lifetime which diverges exponentially as A^c — ^ cxd. 

The region of deconfined plasma metastability terminates when the largest chemical 
potential reaches a maximal value, 

li^UT) = VAT2 + i?-2 X [1 + 0{\)] , (4.26) 

at which point the effective thermal mass of one or more of the scalar fields becomes 
negative. This is analogous to "spinodal decomposition," which determines the limit 
of supercooling or superheating at a typical first order phase transition. However, for 
weakly coupled A/" = 4 SYM, no new equilibrium state exists for yU > /^max- 

5. Summary and Comparison with Dual Gravitational Analysis 

The resulting phase diagram for weakly coupled A/" = 4 supersymmetric Yang-Mills 
theory is sketched on the left side of Figure A confinement /deconfinement transition 
line separates a confining phase at lower temperature or chemical potential from a 
deconfined plasma phase at higher temperature or chemical potential. The location of 
this line is known in the A — > limit but, for the reasons discussed in section 0, the 
order of this phase transition (at small but non-zero coupling) is not currently known. 
If the maximal chemical potential exceeds 1/-R, then no truly stable equilibrium state 
exists. But, at least for sufficiently high temperatures (T ^ 1/-R), the deconfined 
plasma remains metastable beyond ^ = 1/R, with a lifetime which grows exponentially 
as Nc —>■ oo. This metastable region terminates at a spinodal decomposition phase 
boundary when the largest chemical potential reaches a limiting value /imax(^)- The 
A* = fJ'maxiT) boundary line asymptotically rises linearly with temperature with slope 
-\/A, regardless of how many chemical potentials reach /ic(T) simultaneously. 

When T < 1/R, the location of the spinodal decomposition phase boundary is not 
currently known (at non-zero but weak coupling). Neither are order A corrections to the 
location of the confinement /deconfinement transition line. It is possible that the phase 
boundary remains separated from the confinement / deconfinement transition line for all 

^^Note that decay rates via fluctuations involving large excursions of multiple eigenvalues are sup- 
pressed by additional exponentially small factors, relative to the rates of the single eigenvalue processes 
discussed above. As an extreme case, if all eigenvalues make equally large excursions, as in the con- 
figurations shown in Fig. ||c, then the decay rate will be exponentially suppressed by a factor scaling 
as Ar2/A3/2. 
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T T, 

Weak Coupling Strong Coupliing (Gravity) 



Figure 6: The phase structure of the weakly coupled gauge theory on the left, and of gravitational 
solutions to 5-d supergravity, believed to correspond to the strongly coupled gauge theory, on the right. 
For both figures, fi is the largest of the three chemical potentials. In the weak coupling diagram, the 
location of the phase boundary forming the upper limit of the metastable plasma phase has not been 
evaluated at low temperature; this portion of the boundary is indicated by the dashed line. In the 
strong coupling diagram, Tq is the temperature at which the black hole radius reaches to zero. Its 
value depends on the pattern of chemical potentials and equals {ttR)~^ for a single non-zero charge, 
while Tq — for three equal charges. 

non-zero temperatures, with both hnes intersecting at T = and fi = 1/R. This possi- 
bihty is sketched in Fig. ^ But, for non-zero couphngs, the confinement /deconfinement 
transition fine may instead intersect the phase boundary fine at a non-zero tempera- 
ture. In other words, the deconfined plasma region of the phase diagram could pinch 
off before reaching zero temperature. In this case, for sufficiently low temperature, 
the confined phase might extend all the way to the boundary of the phase diagram at 
= ^ma.x{T), with no intervening phase transitions. Which possibility occurs may well 
depend on the particular pattern of chemical potentials (for example, whether two or 
more chemical potentials coincide). 

The right side of Fig. ^ shows a sketch of the behavior of gravitational solutions of 
five dimensional M = 2 gauged supergravity, which are believed to be related to the 
strong coupling limit of A/" = 4 supersymmetric Yang-Mills |]T7| , . The similarities are 
obvious (in contrast to the previous situation depicted in Fig. Q). A qualitatively sim- 
ilar confinement /deconfinement (or Hawking-Page) transition line is present on both 
sides. For the black hole, this transition is known to be first order. The black hole 
instability line, where small perturbations become thermodynamically unstable, ap- 
pears completely analogous to the spinodal decomposition phase boundary of the weak 
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coupling diagram. The black hole instability line rises linearly for temperatures large 
compared to the inverse of the AdSs curvature radius, 1/-R, but the slope depends on 
the pattern of chemical potentials (unlike the case at weak-coupling). With only one 
non-zero potential, for example, the asymptotic slope is vr/y^, while for three equal 
chemical potentials the slope is 2n. 

As noted in footnote ^ of the Introduction, the temperature Tq at which the 
Hawking-Page transition line meets the black hole instability line depends on the pat- 
tern of chemical potentials. For three equal chemical potentials (or equivalently three 
equal charges), Tq = 0. With a single non-zero charge, the temperature Tq = (vri?)"^. 
In every case, the chemical potential at the intersection of the instability and Hawking- 
Page transition lines is 1/R. This intersection corresponds to the point where the black 
hole horizon radius shrinks to zero. Since the phase diagram in the strong coupling 
region is obtained by comparing the values of the action for black hole and AdS solu- 
tions at the same temperature, the region of the phase diagram below Tq is completely 
undetermined. 

Ignoring the low temperature region (where limited knowledge on both sides pre- 
vents a clear comparison), there is one glaring difference between the weak and strong 
coupling phase diagrams of Fig. ^ — the non-perturbative metastability of the decon- 
fined plasma phase when the maximal chemical potential exceeds 1/R. If there is a 
smooth interpolation between weak and strong coupling, then there must be a simi- 
lar non-perturbative black hole instability, with an exponentially suppressed nucleation 
rate whose exponent scales linearly with Nc- No such instability, with an onset precisely 
at fi = 1/R, is currently known in the gravity dual. 



6. Outlook 

There are a variety of possible extensions which should be feasible and which would 
shed light on some of the issues discussed in this paper. 

On the gravitational side, investigation into non-perturbative instabilities of RN- 
AdS black holes (and their generalizations to A/" = 8 gauged supergravity) is clearly 
needed. Our weak coupling results, plus presumed smooth interpolation between weak 
and strong coupling, requires the existence of a non-perturbative black hole instability 
with an onset (at fi = 1/R) prior to the known point of perturbative instability. 

Under gauge/string duality, the eigenvalues of A/" = 4 SYM scalar fields are iden- 
tified with the locations of D3-branes in Type IIB string theory. The i?-charged black 



hole solutions arise from spinning D3-branes [|l^. This suggests that the string theory 



manifestation of the non-perturbative gauge theory instability should involve a stack 
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of D3-branes, for sufficiently high spin, sphtting up into widely separated branes. For 
/i > 1/-R, the force between branes at sufficiently large separation should become repul- 
sive. Since the dominant instability in the gauge theory involves a large excursion of a 
single eigenvalue, the dominant high-spin D3-brane instability should involve a single 
brane separating from the rest of the stack. It should be possible to verify this scenario 
with a probe-brane analysis in the RN-AdS background. 

Alternatively, one may consider the analogue of non-perturbative gauge theory in- 
stabilities involving large excursions by multiple eigenvalues, namely a process in which 
a stack of N^. coincident D3-branes separates into a multi-stack configuration, with each 
stack containing an order one fraction of the N^. branes. (Fig. §c and d correspond to 
the simplest case of a symmetric two-stack instability.) Though subdominant, this 
decay process should be well described by supergravity approximations. 

Finding such an instability for spinning D3-branes, with the expected onset thresh- 
old of /i = would provide further evidence supporting a smooth interpolation be- 
tween weak and strong coupling in A/" = 4 supersymmetric Yang-Mills theory (and the 
validity of AdS/CFT duality). It might also clarify the puzzle, mentioned in Section]^, 
involving the behavior of the gravitational system with large but unequal chemical po- 
tentials, at temperatures below Tq (where the horizon radius of the i?-charged black 
hole shrinks to zero). At the moment, no upper limit on the chemical potential of the 
gravitational phase diagram in this regime in known. 

Turning now to the weak coupling analysis, the order of the confinement /deconfine- 
ment phase transition at non-zero coupling remains uncertain, for the reasons discussed 
in Section 0. Its determination requires a three-loop computation similar to (but more 
complicated than) the one carried out for large Nc pure Yang-Mills theory in Ref. 
This computation, if carried out for non-zero chemical potentials, would also allow one 
to determine how weak coupling corrections shift the location of the transition line, 
which is of particular interest at low temperature where the transition line approaches 
the spinodal decomposition line. If the interpolation between weak and strong coupling 
in AA = 4 SYM is smooth, then the confinement /deconfinement transition is expected 
to be first order. Confirming (or refuting) this directly would clearly be desirable. 

Finally, our weak coupling analysis focused on the high temperature regime, T 3> 
1/ R. Extending the treatment to T < 1/ R would be valuable. For T <^1/ R one has a 
large time circle and a small spatial sphere, so the non-zero angular momentum modes 
on the three-sphere are heavy and may be integrated out, leading to an effective one 
dimensional quantum mechanics. Analysis of the resulting effective theory should reveal 
the behavior of the spinodal decomposition phase boundary in the low temperature 
region, determine whether it intersects the confinement /deconfinement transition line 
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at a non-zero temperature Tq, and perhaps even reveal a distinct high density, low 
temperature phase. 
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A. Derivation of Matrix Model with Chemical Potentials 



To reduce the partition function of the free theory (constrained by Gauss' law) to the 



matrix model, as sketched in section |2.3| , one first separates the zero mode of Aq by 
replacing 

A,y{x) — > A,y{x) + 6„fl a/g , (A.l) 

where a is a constant traceless Hermitian Nc x Nc matrix, and Aq is orthogonal to 
the constant mode. Without loss of generality, we may assume that a is diagonal and 
we denote its entries (or eigenvalues) by {q"^\m = 1, ■ ■ ■ ,Nc} where the g™" are all 
real and satisfy = 0. One plugs this shifted form of the gauge field into the 

Lagrangian ( |2.19| ) and then sends g —>■ 0. It is convenient to choose the gauge defined 
by the gauge-fixing term, 

tT{{d,A, + t[a,Ao]f}. (A.2) 

Define, for convenience, g™" = q"^ — q"- and D'"" = {di, + iq^^5pflY . The Lagrangian 
of the free theory may be rewritten as the quadratic form 

-^quad = {■^y)*mn{~^^^){-^i')mn 
m,n ^ 

~^ ((^p)mn? 0^p)mn) 

+ (V'i):„n(#- Tog"'" - ^^^aol5){^^)mn } , (A.3) 



where the color indices m,n are now explicit, with {Ay)mn = A^ (T°')mn and likewise 
for the other fields. 

Table |^ in section |2.3| lists the eigenvalues of the spatial small fluctuation operators 
whose eigenfunctions are (scalar, vector, or spinor) spherical harmonics. Denote these 
eigenvalues by A^, and A^ for the transverse spatial gauge, conformal scalar, and 
spinor fields, respectively. Denote the Matsubara frequencies by ook, for either bosons 
or fermions (as determined by context). Then the contribution to InZ from the gauge 
bosons and ghosts (equal to — | times the logarithm of their functional determinant) is 



5^ 5^ trln{(c.. + q-r + A^} - tr ln(u;^ + Aj) 



>- k 



(A.4) 



The trace is over the space of transverse vector fields on S^. The contribution from 
gauge-fixing ghosts cancels the contributions from Aq and longitudinal part of the 
spatial gauge field, thereby effectively leaving two transverse gauge field degrees of 
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freedom. (This factor of two is included in the degeneracy factor for transverse vector 
fields listed in Table |[) The second term compensates for the fact that there are really 
only Nc—l independent diagonal components of the gauge field, not N^, because we 
are dealing with the gauge group SU{Nc). 

Block-diagonalizing the linear operator in ( |A.3|) for the scalars (by working in the 
eigenspaces of do and the spatial Laplacian) and taking the determinant in each block 
gives the scalar contribution of 



k p '- m,n 



— tr In 

Similarly, the fermions give 



(A.5) 



k i 



m,n 



tr In 



{{ujk + ii^iY + A;}{(cufc - ifii)' + A^} 



(A.6) 



In integrating out the fermions, it is essential to recall that ipi is a Majorana fermion. 
So ipi and ipi are not independent and the Grassmann integral gives the Pfaffian of 
(^+ i7o g™*^ — /iiToTs), not the determinant. 

The Matsubara frequency sums can be carried out (most easily by differentiating 
with respect to A^, performing the sum, and then integrating back in A^) and yield 

InZg = -i ^trln[{l - e-'3(A«+*'?'"")||i _ e-^^^^^"™")}] + tr ln(l - e'^^^ , (A.7a) 



InZs = - ^ 1 5^tr ln[{l - e-^(^=+^-+^'^'"")}{l - e-/5(^^+^.-'?-)} 



m,n 



X {1 - g-/3{As-Mp+ig™")j.|;L _ g-/3{A.-Mp-ig"")| 



tr In 



(A.7b) 



i '- m,n 

X {1 + g-/3(A/-M»+ig™")}|l _^ g-/3(A^-M«-ig™) j 

- trln[{l + e-^(^/+^>)}{l + e-/3(^/-'^')} 



(A.7c) 
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In these expressions, we have dropped temperature independent terms that do not 
contribute to the thermodynamics of the theory. To evaluate the remaining traces 
(over transverse vector, scalar, or spinor fields on 5*^), it is convenient to expand each 
logarithm in a power series. For example (taking part of the scalar contribution) 

-Yl iE^^Ml - e-'^(^=+^''+^'^'"")} - itrln{l - e-^(^^+'^-)} 

p L m,n 

oo oo ^ / 

p h=0 1=1 \m,n 

1=1 \h=0 / \ p / \ m n ) 

= \ t (T. U' C/t' - 1) , (A.8) 



where x = e ^ = e ^^"^ and U = exp[i/3a] is an Nc x Nc unitary matrix. Similarly, one 
can work out the other logarithms. One finds 



In = y ^ ^ =^ (tr tr t/t' - l) 

1=1 ^ ' 



, ^ 1 x' + x^' 

In Zo = > - — ■ 

/ (1 



,p=l 



oo 



/ {1-x^y^ 

/ (l-x^\^ 



J2 i^^^' + 



X 



.1=1 

■ 3 

n 

.p=i 



(A.9) 

(tr t/' tr [/^' - 1) , (A. 10) 

(tr f/' tr f/"f' - 1) 

(tr f/' tr f/^' - 1) , (A. 11) 



where the definition ( |2.12| ) of the effective fermion chemical potentials /Xj was used in 
the last line. 

Combining these terms and defining 



5'efr(f/) = — In — In — In Z 



(A.12) 



gives the result quoted earlier in (|2.25|) and ( p.26|) . The exponential e gives the 

result of integrating out all fields except the zero mode of A^. What remains is a single 



integral over the matrix [/, as shown in (|2.27|) . 
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B. One-Loop Matching for ESYM3 



B.l CoefRcient of the Identity Operator 

The leading contribution to the coefficient of the identity operator in ESYM3 may 
be obtained by computing the functional determinants which resuh from integrating 
out the heavy modes, neglecting interactions. To evaluate the contributions from the 
heavy bosonic modes, ffist consider an 0(2) invariant theory (in flat space) with two 
real massive scalar fields. One may easily evaluate the contribution to the free energy 
produced by integrating out these fields. With a non-zero chemical potential, the high 
temperature expansion of the result is 



F ^2_2^2 ^2 2 3/2 

- = + t^T^-S tl2 — r + --- , B.l 

V 45 12 671 ' ^ ' 

where V is the volume of the space. This result contains the contributions of both 
the heavy non-zero Matsubara modes and the light zero-frequency mode. For our 
purposes, the light mode contribution must be subtracted out. The zero Matsubara 
mode contribution to the above expression is given by the integral, 

- T —\:r ln{p^ + - fi^) = TT-^ — (B.2) 



V J (27r)3 ' Qtt 

This is exactly the third term in F/V. This was inevitable — such a non-analytic term 
is related to the infrared behavior of the theory and can only come from fiuctuations 
on scales large compared to 1/T. Therefore, the heavy mode contribution is 



V [45 12 V ^2 J \^ ' ' y4 ^ 

In A/' = 4 SYM, we have 3N^ such pairs of scalars. Moreover, gauge bosons and ghosts 
contribute times the above expression with m = yU = 0. Hence the total bosonic 
contribution is. 



_ ^2^4 I ^ 



V ' 45 12 



0{\') , (B.4) 



where we have set iv? = 1/ R^. Because we are working in the regime where l/{TKf' ~ 
Ijp' jT"^ ~ A, this result is correct through 0{X^). 

The analogous contribution from one Majorana fermion with chemical potential /i 
[introduced as in Eq. ( |2.14| )] is 



V V360 12 r2 24 ' ^ ' 
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Since fermions do not have zero-frequency modes, there is no hght mode contribution 
to subtract. We have altogether AN"^ Majorana fermions with four different chemical 
potentials, so the total one-loop fermion contribution is 

= -N^X + -L { i^^^^^i^ ^ o(A-)} , (B.6) 

Combining the boson and fermion results, we obtain the one-loop contribution to 
the free energy density, expanded in powers of \/{TR) and yu/T through (9(1/T^), 

^i-ioop _ NlTU . 3 , 2{^^\ + ^^l + ^^l) , + + fil + fil , 

-^--^Y"" "r^ + + + 

(B.7) 

However, this is not a complete result. Because we are interested in the regime where 
\/{TRY ~ ^J-^/T'^ ~ A, two- loop contributions will yield C(A) corrections to the free 
energy density which are comparable in size to the curvature and chemical potential 
dependent parts of the one-loop result. The required two-loop contribution may be 
evaluated ignoring curvature and chemical potential corrections altogether, since these 
will only give subleading 0{\/ {TRY) or (9(A/i^/T^) terms, which are both order A^. 
One finds [||, 11,11, 



= ^ NlT^ A + 0{\^) . (B.8) 

With massless fields in fiat space, there are no light-mode contributions which need to 
be subtracted out. Thus, the coefficient of the identity operator in ESYM3, equal to /3 
times the free energy density due to the heavy modes, is given by 



/ = -^ 2^^-3A-^+ "--^ ' -^- + -^ -^ +0(A 



B.2 Scalar Thermal Mass 



(B.9) 



The lowest-order thermal mass correction for the scalar fields comes from the self-energy 
diagrams shown in Fig. |^. As noted in section ^ to extract the (9(AT^) contribution to 
the mass (squared), one may evaluate these diagrams in fiat space and without chemical 
potentials. Including curvature and chemical potential corrections in these diagrams 
would give subleading O(A^) terms, which are beyond the accuracy we need. Since 
we are interested in the construction of ESYM3, only the heavy mode contributions to 
each loop should be included. To obtain the thermal mass, it is sufficient to evaluate 
these diagrams at zero external momentum. But we will keep the external (Euclidean) 
momentum p = (0,p) non-zero (but small compared to T) in order to extract the 
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(51) 



(S2) 



(S3) 



(54) 



(S5) 



Figure 7: One-loop self-energy diagrams for the scalar fields. The solid lines are scalars, curly 
lines are gauge bosons, and dotted lines are Majorana fermions. The final diagram (5*5) denotes 
the wavefunction renormalization counterterm. The external lines are light modes with small spatial 
momenta, while the fields running inside the loops are heavy modes. 

one-loop wavefunction renormalization for the scalars, which will be useful in the next 
subsection. We use dimensional continuation and the MS renormalization scheme, 
and denote the arbitrary renormalization scale by A. Integrals are carried out in the 
standard manner with the convention that 



The prime in the left-hand side signifies the omission of the zero frequency mode in 
bosonic loop integrals. We use a Lorentz gauge-fixing term with arbitrary gauge pa- 
rameter a. The dependence on a will drop out of the result for the thermal mass 
correction, but keeping a arbitrary serves as a useful check on the calculation. 




(B.IO) 



Evaluating the diagrams of Fig. |^ is straightforward and one finds: 





(B.llb) 




(B.llc) 



(^5) 





45 



where SZ^^^ is the one-loop wavefunction renormahzation counterterm. We have defined 



Lb = ln^-21n(47r) + 27E, = In ^ - 2 ln(47r) + 27^ + 4 In 2 , (B.12) 



with jE is Euler's constant. The sum of these diagrams gives (minus) the hard mode 
contribution to the one-loop self-energy, 



-nhard(p) = iSl) + iS2) + {S3) + (SA) + (55) 

= -^T^+p"y^2 {(3-«)L, -4L, - i±^} + 0(| 



4 



(B.13) 



Choosing 

5Z« = -(l + a)^ (B.14) 

absorbs the logarithmic divergence. As always, temperature independent renormahza- 
tion suffices to remove UV divergences at any temperature. The hard-mode self-energy 
is thus 

-nhard(p) = -XT' + p2_A_ {{3-a)L, - ALf} + o(^^). (B.15) 

The resulting renormalized scalar propagator, in fiat space and without chemical 
potentials, at separations large compared to 1/T, is 



{^Aix)^Biy)) = ^abT j e^-^-^ Gip) , (B.16) 



with 



G-\p) =p' + nhard(p) + nsoft(p) 

= p\l + A) + XT' + nsoft(p) + Oip'/T') + 0{X'T') , (B.17) 

where A = {{3—a)Lh — ALf} and nsoft(p) is the self-energy contribution due to 

soft modes. 

In view of the relation ( p.l| ), the long distance correlator ( p.l6|) should be matched 
with Z3T($a(x)$b(|/)), computed in the effective theory. Again neglecting curvature 
and chemical potential corrections, the effective theory correlator is 

d^p 



^A{x)^B{y)) = Sab 1 ^ e^'^-^'-y^ G{p) , (B.18) 



with 

G-\p) = p2 + Srn'(^T) + U^^M , (B.19) 
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and nsoft(p) is the same soft mode contribution as in (|B.17|) . Therefore, the parameters 
of the effective theory must be adjusted to make G^^{j))/Z2, coincide with G~^{j)). 
Through one-loop order, the required matching is 



A 



1 + A)-' = 1 + {(3-«)L, - 4L^} + 0{\') 



and 



6m\T) = XT' + 0{X'T'). 



(B.20) 



(B.21) 



B.3 Scalar Quartic Coupling 

Calculating the one-loop thermal corrections to the scalar quartic interactions which 
will appear in the high temperature effective theory requires evaluation of the one- 
loop four-point correlator of light-mode scalar fields, at zero external momentum. The 



analysis is similar to the computation by Nadkarni of the quartic interactions of Aq 
in EQCD. The required diagrams are shown in Fig. ^. Once again, only heavy modes 
should be regarded as running around the loops, and one may work directly in flat 
space, ignoring curvature and chemical potential corrections. 




Figure 8: One loop diagrams which lead to quartic scalar interactions in the high temperature 
effective theory. As in Fig. solid, curly, and dotted lines denote scalar, gauge boson, and fermion 
lines, respectively. In Diagram {Q6), it is most convenient to treat the external fields X and Y 
separately due to their different Yukawa couplings. Diagram (Q?) denotes the contribution induced 
by wavefunction renormalization. Not shown explicitly are the various permutations of each diagram. 
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We evaluate these diagrams using dimensional continuation and MS renormaliza- 
tion, as described in appendix p.2| , and an arbitrary gauge fixing parameter a, even 
though choosing Landau gauge {a = 0) would reduce the number of the diagrams to be 
computed (namely diagrams {Q3), ((54), and (QS) would all vanish at zero external 
momentum). Nevertheless, retaining an arbitrary gauge fixing parameter in order to 
verify the cancellation of a dependence in the resulting effective quartic coupling serves 
as a useful consistency check. 

To express the results compactly, we will use {V^'^} to denote Hermitian SU{Nc) 
adjoint representation basis matrices, with^^ 



(V- 



be 



abc 



(B.22) 



and define the following index structures: 



fj^ \abcd 
k'- IJABCD 



n: 



T \ abed 
2JABCD 



= n: 



triV'^V'VV'') i5AB ScD + Sad Sbc - 25 ac Sbd) 
+ ti{V''V''V''V') {6 AC Sbd + Sab Scd - 26 ad Sbc) 
+ tiiV^V'V'V'') {6ac Sbd + Sad Sbc - 26 ab 6cd) 

tiiy'V'V'V'') {5ab Scd + Sad Sbc) 
+ tiiV^V'V'V') {6 AC Sbd + Sab Scd) 
+ tiiVV'V'V'') (Sac Sbd + Sad Sbc) ■ 



fB.23a) 



;B.23b) 



The tree-level quartic scalar vertex is —2 A (Ti)^'']^'^^ . 

The diagrams of Figure H can be computed in a straightforward manner, and one 



^^Useful adjoint representation trace identities include 

tr(y"y^) = Nc 5"^ , 
tr(y"T/V) = -|^^c^''^ 

tr(F"y V^F'') = Nc tr{T''T''T''T'^ + t^t'^T'^T') + \ {S'^^S"'^ + 5'"'5'"^ + 5'"^5'"' 
where {T°} are the fundamental representation basis matrices. 



(B.24a) 
(B.24b) 
(B.24c) 
(B.24d) 
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finds 



(Ql) = (7 + Lt) [m)ABCD + mnBco] , (B.25a) 

(Q2) = (3 + a') + L,) iT,)%'cD , (B.25b) 

(Q3) = -2a2 (7 + (r2)t^'cD , (B.25c) 

(Q4) = (7 + ^^.) (r2)t^'cD , (B.25d) 

m = -Aa Y^ (7 + U) {Tr)%'cD , (B.25e) 

(Q6) = -8 Y^ (7 + ^/) [(ri)^^'cD + (r2)^|'cD] , (B.25f) 

(Q7) = -4A5Z«(ri)^^"cz.- (B.25g) 



Diagram (Q6) requires a special attention because the Xp and 1^ scalars have different 
Yukawa couphngs. Considering the various cases (such as {X^ X'^Y^Y^)) and taking 
care of the appropriate permutations, one finds the simple expression above. Inserting 
the previously determined one-loop scalar wavefunction renormalization factor 5Z^^^ 
[from Eq. (iBlil) ] into the result (Q7) gives 

m = 4 (1 + a) {T,)%i,^ . (B.26) 

Combining all terms, we find that the amputated IPI scalar 4-point function, at 
zero momentum, is 



+ {^2^300 [Lt - Lf] + 0{X') . (B.27) 

Rescaling the field according to relations ( |3.1| ) and (|B.20| ) removes the remaining gauge 
fixing parameter dependence and produces the equivalent ESYM3 IPI 4-point function, 

(^^^b^c*d)ipi = -2AT2 ir^)%%^ - (In 4) ^ [{T2)%'co " i^i)ABCD] + 0{\') • 

(B.28) 

The first term is the vertex which arises from the commutator squared interaction 
directly inherited from the full theory, while the order correction comes out to be 
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remarkably simple. Both terms must be produced as tree level vertices in the three- 
dimensional effective theory. The required quartic interaction terms in the ESYM3 
effective Lagrangian are 

\gl tr(^[|.^, <l^])2 + ^ I iA^WAy^f , (B.29) 

with g\ = g^T . All gauge fixing parameter dependence and 1/e poles have canceled, 
as required. The ln(A^/T^) terms in Lf, and Lf have also dropped out. This was 
also required, since A/" = 4 SYM at zero temperature is a conformal theory in which 
the coupling does not run, and turning on a temperature does not change the short 
distance behavior. 



C. Background Field Method 

We wish to evaluate the complete one-loop effective potential for the scalar fields using 
background field techniques. We will restrict our analysis to the special case where 
the background scalar fields take values along fiat directions of the tree-level potential. 
This is sufficient for our purposes, and simplifies the calculation. We carry out the 
calculations in flat space, and neglect chemical potentials. Under our assumptions that 
1/ and /i^ are both of order AT^, the tree- level contributions from the curvature and 
chemical potential induced mass terms will be comparable to one-loop corrections to 
the effective potential computed in flat space with /i = 0. Including curvature and 
chemical potential corrections in the evaluation of the one-loop potential would change 
the result by an amount comparable to neglected two-loop contributions. 

In the flat space Lagrangian, shown in Eq. ( 2.19| ) but without mass and chemical 



potential terms, we shift the scalar fields by constants, 

^A-^^A + <yA , (C.l) 

with $^ constant fields taking values in the flat directions, and cTyi, are arbitrary fluc- 
tuations about the constant background. The background fields (viewed as Nc x N^. 
matrices) are simultaneously diagonal (up to an irrelevant gauge transformation). 

Among the quadratic terms in the shifted Lagrangian, there are terms which mix 
the gauge and scalar fields. These may be eliminated by employing the -R^-gauge with 
gauge-fixing term 

-^ii{{d,A^-iig[aA,^A]?] . (C.2) 

We further set ^ = 1 to simplify the computation, as this choice eliminates off-diagonal 
quadratic terms which mix aA and (Tb for A ^ B. We denote the matrix components 
of the fields by indices m and n so that, for example, (o"^)mn = Sa'^l(^'^)mn- 
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Let denote the mth eigenvalue of the constant background field $a and define, 
for convenience, 

Ml^^g'Y.^^A-K,f ■ (C.3) 



A 



Note that m and n are not to be summed over in the definition of M^„. The resulting 
quadratic terms in the shifted Lagrangian are 



(C.4) 



In carrying out the Gaussian integrals, one finds that the logarithms of the functional 
determinants from all the three contributions simplify to the generic form 



ln(p^ + M^J, (C.5) 

where the sum-integral is defined in Eq. ( [B.1CI|) . This integral is standard after Arnold 
and Espinosa [^, and yields 

^22^4- 45 + 3[^^2T2j VUTT^T^y (2/)!! '\An^T^ 



(C.6a) 



^2rp4 360 e\4T,2rp2j 2^^' ^JUn^TV ' (2/)!! 'K^n'^T^ 

(C.6b) 

where /bose and /fcrmi are contributions from the bosons and fermions, respectively. In 
these results, 0(e) terms have been discarded and the quantities and Lf are defined 
in Eq. (|EA^ ). 

For each color component, the gauge field, six scalar fields, and gauge fixing ghosts 
contribute (4 + 6 — 2) = 8 times f/bose to the (three dimensional) effective potential, 
while the fermions contribute —8 times f /fermi- The resulting effective potential, for 
background fields lying along flat directions, is thus 

F(S,) = l.=r' 5: [ - 1 + (^) + (In 2) (^|§) ' 



2 \ ' 

■ • (C.7) 
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The 1/e poles and all dependence on the arbitrary renormalization scale has canceled, 
as required for a physical quantity. 

We would like to express the result ( |C7] ) in terms of the original fields $a- To do 
so, note that the adjoint representation of the field, {^^^)ab = may also be 

written as 

($f U = 2tr(Tj$^,T,]). (C.8) 

Since the fields $a are assumed to take values in flat directions, a Cartan subalgebra 
may chosen for which $^ is diagonal with the entries A^. Therefore 

a 

= 2 ^ ^ (^a)nm(A™ — A^) (A^ — A^) (TQ,)mn 
a,m,n 

= E(^^-^^)(Ab-AS), (C.9) 

m,n 

where the last step used the identity J2ai'^o.)ij{Ta)ki = liSuSjk - -^^kjhi)- Hence 

W Ea ^f^f] = 9' EA,m,niK-Kr = Em,n ^L- Siuce the sum of the eigenval- 
ues A™ must vanish (because $^ is traceless), this may also be written in the form 

Ig'T'tri^f^f) = iT^J^ML =/m$^(A:?)^ = XTHri^A)' , (C.IO) 

m,n A,m 

which shows that the quadratic part of the effective potential ( |(J.7[ ) agrees with the 
previously derived thermal mass correction ( p.6|) . [Note that, to the lowest order, the 
rescaling from $^ to $a in (|3.1| ) is just a factor of VT.] More generally, we have 



tr 



A m,n 



SO the effective potential ( |C.7| ) (along flat directions) may be expressed in terms of 
adjoint representation traces as shown in Eq. ( |3.11|) . In particular, the quartic term 
of ( |(J.7| ) is (In 2) (jf^ tr[(<l>^^-'<l>^^-')^] and this agrees, for commuting fields, with the 
diagrammatic result ( [B.29| ).^^ 
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The meticulous reader may be puzzled why this result agrees with Eq. (B.IO), without taking 



into account the order A term in the wavefunction renormalization Z3 [given in Eq. (B.20)]. As seen 



explicitly in Appendix B.3, this rescaling leads to a term proportional to the tree- level interaction and 



hence vanishes for fields lying along flat directions. 
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